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Abstract 

We use the consistency approach to classify discrete integrable 3D equations of 
the octahedron type. They are naturally treated on the root lattice Q(A^) and are 
consistent on the multidimensional lattice Q(An). Our list includes the most promi- 
nent representatives of this class, the discrete KP equation and its Schwarzian (multi- 
ratio) version, as well as three further equations. The combinatorics and geometry 
of the octahedron type equations are explained. In particular the consistency on the 
4-dimensional Delaunay cells has its origin in the classical Desargues theorem of pro- 
jective geometry. The main technical tool used for the classification is the so called 
tripodal form of the octahedron type equations. 
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1 Introduction 

This paper is devoted to the study of certain integrable discrete three-dimensional equa- 
tions. Prototypic examples of such equations are the discrete KP (dKP), or the discrete 
bilinear Hirota equation, and its Schwarzian version. Both equations play a pivotal role in 
modern mathematics [16, 18, 21, 22, 23, 24, 25, 26, 27, 31, 36, 39]. 

The dKP equation for a complex lattice function (field) x : Z 3 — > C was introduced in 
[20] in the form 

ax(m + ei)x(m — e±) + f3x(m + e^xiyn — e 2 ) + 72(771 + e 3 )x(m — e 3 ) = (1) 

(here m £ Z 3 , ej is the unit vector of the i-th coordinate direction, and a, f3, 7 are three 
arbitrary complex coefficients). It relates six fields assigned to any elementary octahedron 
of the even sublattice Z 3 ven of Z 3 , centered at the odd point m £ Z 3 . The even sublattice 
Z 3 vcn is also known as the face centered cubic lattice (fee lattice). Its 3-cells are octahedra 
and tetrahedra. 

The fee lattice is isomorphic to the root lattice of the type A3, 

Q(A 3 ) = {n= (n , . . . , n 3 ) G Z 4 : n + m + n 2 + n 3 = 0}. (2) 

The six vertices of every octahedron cell in Q(A 3 ) can be written as n + e% + ej, i,j £ 
{0, 1, 2, 3} (with n from the hyperplane uq + n\ + n 2 + n 3 = —2). The dKP equation relates 
the fields x at these vertices and reads: 

ax(n + e + e\)x(n + e 2 + e 3 ) + f3x(n + e + e 2 )x(n + &\ + e 3 ) 

+7x(n + e + e 3 )x(n + e\ + e 2 ) = 0. (3) 

This representation of the dKP equation turns out to be more conceptual and plays the 
key role in the integrability analysis of this paper. 
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Since on the hyperplane uq + n\ + n 2 + ra 3 = any coordinate is determined by 
other three, by forgetting one coordinate we obtain a one-to-one correspondence of the 
lattice Q{A$) with Z 3 , (no, n\, n 2 , n 3 ) (ni, rt2, «3)- In this realization, the dKP equation 
becomes an equation on Z 3 : 

ra(n + ei)x(n + e 2 + e 3 ) + /3x(n + e 2 )z(ra + e 3 + ei) + 7x(n + e 3 )x(n + e\ + e 2 ) = 0. (4) 

In this form, introduced in [28], it relates fields assigned to six out of eight vertices of any 
elementary cube of Z 3 (the fields at one pair of the opposite vertices, n and n + e\ + e 2 + e 3 , 
do not participate in the equation). 

Throughout the present paper, we will use the following abbreviation for lattice func- 
tions: x for x(n), Xi for x(n + ej), x-i for x(n — e^), Xij for x(n + ej + e^), etc. In this 
notation, the dKP equation in the Q(^4 3 )-form (3) can be put as 

aX lX 2 3 + /3^02^13 + 7^03^12 = 0, (5) 

while in the Z 3 -form (4) it can be put as 

axix 2 3 + /3x 2 x 13 + 7x 3 x i2 = 0. (6) 

We will call equations like (6) the octahedron type equations, as opposed to the cube type 
equations, whose best known representative is the dBKP equation introduced in [28]: 

UX1X23 + /3x 2 xi3 + 7X 3 Xi 2 + 5XX\23 = 0. (7) 

This latter equation relates the fields at all eight vertices of any elementary cube of Z 3 . The 
octahedron type equations could be considered as a subclass of the cube type equations, 
however they have different properties and require for a different analysis. 

Definition 1. Equation of octahedron type is the relation 

F{x i, x 02 , x 03 , x 12 , x 13 , X23) = 

for the unknown function x : Q(A 3 ) — > C, or, equivalently, the relation 

F(xi,X 2 ,X3,X 12 ,Xi3, X23) = 

for the unknown function x : Z 3 — > C. 

Our only assumption concerning the function F is that it should be locally analytic 
and satisfy 

irreducibility condition: equation F = can be locally solved with respect to 
any variable, the result depending on all other variables, i.e., on any solution 
of F = we have F x ^ , where x denotes any argument of the equation. 

This condition forbids, in particular, equations with F of the form F = AB, where A 
or B does not depend on some variable (for example, if A x = 0, then F x = AB X , and 
this vanishes on the solution A = 0). Also the so called ultradiscrete equations with 
piecewise constant functions F are excluded from consideration. On the other hand, the 
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irreducibility condition does not forbid the case when the solution of F = is multivalued, 
for instance, F = AB with both A and B depending on all variables. We handle such 
situations by working always in the neighborhood of some solution branch, where the 
theorem on implicit function applies. 

It is important to observe that irreducible equations of the octahedron type are 3D 
systems in the following sense: a generic solution of such an equation on Z 3 can be defined 
by a well posed Cauchy problem with generic 2D initial data. For instance, for the dKP 
equation in the form (1) such initial data are constituted by prescribing the values of x on 
the planes n% = and n% = —1. 

In the work [3] we pushed forward the idea that integrability of discrete equations is 
synonymous with their multidimensional consistency. For a certain class of 2D equations, 
this notion was put in the basis of a classification of integrable cases. In the present work, 
we classify multidimensionally consistent (integrable) equations of the octahedron type. A 
general idea of consistency leads to the following formulation. 

Definition 2. Consider a system on Z N consisting of (possibly different) octahedron type 
equations 

F(xi, Xj , Xj~i Xij , Xik 7 Xjfc^) — (8) 

on all affine 3D sublattices c+Zei+Zej+Zek- It is called multidimensionally consistent if it 
has a solution whose restrictions on all 3D sublattices are generic solutions of corresponding 
equations. 

This definition is a literal repetition of the corresponding notion for cubic type equations 
and does not take into account the specific feature of the octahedron type situation. If the 
lattice Z N lattice in this definition is treated as a realization of the root lattice 

Q(A N ) = {n = (n , . . . , n N ) G Z N+l : n + m + . . . + n N = 0}, (9) 

by forgetting the coordinate no, i.e., Q(A^) 3 (no, m, . . . , n/v) (ni, ■ ■ ■ , njv) G Z N , then 
the corresponding octahedron type equations (8) live on the sublattices 

Q(A 3 ) = {(no, rii, nj,n k ) : n + n, + nj + n k = const}. 

However, these are by far not all Q(A^) sublattices of Q(A^), but only those involving the 
coordinate no, which becomes distinguished in this formulation. A more symmetric notion 
would be: 

Definition 3. Consider a system on Q(An) consisting of (possibly different) octahedron 
type equations 

on all affine sublattices Q(A^) = {(ni,nj,nk,n m ) : m + rij + n k + n m = const}. It is 
called multidimensionally consistent if it has a solution whose restrictions on all Q(A%) 
sublattices are generic solutions of corresponding equations. 

This version is more natural but looks much more restrictive, since there are ( N ^ 1 ) 
sublattices Q(A%) of Q(A^), and only ( 3 ) sublattices Z 3 of Z N . For instance, in the 
case N = A the lattice Q(A±) has 5 sublattices Q(As), of which only four involve the 
distinguished coordinate axis uq. Nevertheless, in Section 3.5 we prove the following result 
(see Proposition 4). 
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Theorem 1. Definitions 2 and 3 of the multidimensional consistency of octahedron type 
equations are equivalent. 



Now we formulate the main result of this work. We classify multidimensionally consis- 
tent systems of octahedron type equations (in notation of Definition 3) modulo the group 
of admissible transformations. This group consists of changes of independent variables n 
generated by the affine Weyl group of Q (A jy) (permutations of indices and translations by 
lattice vectors) extended by the simultaneous inversion of all coordinates n i— > —n, as well 
as of non-autonomous point transformations of dependent variables, x(n) h-> f(x(n),n). 

Theorem 2. Any multidimensionally consistent system of octahedron type equations on 
Q{An) can be reduced by an admissible transformation to a system whose restriction to 
any sublattice Q(A^) has one of the following forms: 

XijXkm X-ifcXjni + XjkXi m — 0, (Xl) 

{xirn Xij^(Xj m 2 ; jfc)(^A;m Xik) ^ ^ ^ 

{Xij X jm)(,X jk Xkm)(,Xik Xi m ) 

(Xik Xij)Xi m + \X%j Xj^Xjm + (xjk Xij t )xj {m — 0, (X3j 

Xik X{j ^ Xij Xjk ^ Xjfc Xik q ^ ^ 

Xim Xj m X\ tm 

Xik Xjk / 1 1 \ . . 

— =Xij{ )• (X5) 

Xkm ^ Xj m Xi m / 

Moreover, each of equations (xi), (X2) is multidimensionally consistent with itself (that is, 
with the set of the same equations for all Q(A^) sublattices), while multidimensional sets 
involving any of the equations (X3), (xa) include with necessity (x%) on some of the Q{A%) 
sublattices, and multidimensional sets involving (xh) include with necessity (xa) on some 
of the Q(Az) sublattices. The detailed description of consistent sets is given in Theorem 
21. 

Remark 1. All these equations in Z 3 -form already appeared in the literature. They were 
derived by the direct linearization method in [9, 32, 31, 30, 14, 8]. First steps towards 
consistency of the dKP equation on a root lattice of type A were made in [35, 34]. 



Remark 2. It is instructive to look at the (somewhat more traditional) Z 3 -form of equations 
from Theorem 2. Recall that they are obtained by forgetting one of the indices. The high 
symmetry grade of equations (xi), (X2) yields that forgetting any one of the indices leads 
to the same equation on Z 3 , namely 

XijXfc XifcXj -\- XjfcXi 0, 



resp. 

{Xi Xij){Xj XjkjiXfo 2^jfc) 
(X{j Xj)(Xjfc 2-fc)(^ifc Xi) 
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As for equations (X3)~(X5)> one gets for each of them several seemingly different equations 
on Z 3 , like 

%ik %ij , Xij %jk Xjk %ik q 

Xi Xj Xfc 

and 

H — H = 0, 

Xfc Xjk Xik 

both of which follow from (x<i) by forgetting one of the indices. Without the unifying and 
symmetric Q(A^) notation, it is not easy to recognize the equivalence of the latter two 
equations. For instance, the (non-commutative versions of) these equations are listed in 
[30] as different equations (1.6) and (1.9). 

Remark 3. One can get equations (x-l)) (xs) from (X2) via simple limiting transitions. Per- 
forming in (X2) a non-autonomous point transformation x(n) *— > e~ nm x(n), which amounts 
to the replacement of Xi m: Xj m ,Xk m by e -1 Xj m , e _1 Xj m , e~ 1 Xk m , and then sending e — > 0, 
we arrive at equation (x±j- Analogously, performing in (X4) the point transformation 
x(n) 1 y 5 nk x(n), which amounts to the replacement of Xik,Xjk,Xk m by 8xi m ,5xj m ,8xkm, 
and then sending S — > 0, we arrive at equation (xs)- Of course, these limiting transitions 
do not belong to our group of admissible transformations. 

The contents of the paper is as follows. 

In Section 2 we give, for the sake of completeness, a simple and well known derivation 
of the dKP equation and the related octahedron type equations from the compatibility 
of auxiliary linear problems. The main result of the present work, Theorem 2 (or its 
detailed version Theorem 21), says that equations derived in this section exhaust the list 
of multidimensionally consistent equations of this type. 

Section 3 is devoted to the definition of multidimensional consistency for the octahedron 
type equations. This definition is not quite straightforward, since the underlying lattice 
Q(As) is more sophisticated than the standard cubic one. The main problem is to find a 
suitable multidimensional lattice containing several copies of the Q(A^) lattice which can 
simultaneously support generic solutions of the discrete octahedron type equations. The 
fact that this problem is not trivial is illustrated in Section 3.2 by a failure of one possible 
definition. A successful definition is then illustrated by the example of the dKP equation in 
Section 3.3 and formulated in full generality in Section 3.5. An elementary combinatorial 
cell (a 4-cell of the root lattice Q{A±) with five octahedral faces) is best illustrated by the 
Desargues configuration, see Section 3.4. 

Section 4 contains the central technical observation: each octahedron equation of the 
multidimensionally consistent system can be written in eight ways in the so called tripodal 
form, and the tripodal forms of equations on the adjacent octahedra must combine them- 
selves in a very special way. Actually, the necessary conditions for consistency established 
in this section are the basis for the subsequent solution of the classification problem. 

This solution starts in Section 5 where we classify all octahedron type equations ad- 
mitting eight tripodal forms. This conditions turns out to be stringent enough to produce 
a finite list of equations, given in Theorem 7. 

Finally, in Section 6 we combine tripodal equations into consistent systems, a complete 
list of which is given in Theorem 21. 
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2 Equations of octahedron type through linear problems 



All equations from the list are related to each other via difference substitutions, and can 
be easily derived from simple linear problems like 



h - fi = a /> h~ h= b f- 



(11) 



It should be noticed that all four faces of the tetrahedron with the vertices /, /i , / 2 , /3 are 
on the same footing, because, due to (11), there hold also the further linear equations: 



h ~ h = cf, h - h = d{h - fi 



(12) 



where c = b — a, d = a/b. Any two out of the four equations (11), (12) are equivalent to 
(11). 



h 





i /l3 










b \ 


h \ 




J \ 


fi 


\/n 




ai 











Figure 1. Linear equations on triangles. 

It is easy to compute the compatibility conditions (which express the equality of two 
alternative expressions for /23 through the initial values /, /i, /n, see Figure 1): 



03 + 61 = 62 + ^1, a 3 6 = 6 2 a. 



(13) 



These formulas should be interpreted as a map (a, ai,b, 61) h-» (03,62)- This map defines 
the evolution of the initial data consisting, for instance, of the values of a prescribed on the 
coordinate plane 12 and the values of 6 prescribed on the coordinate plane 13. Expressing 
in (13) a and 6 through / according to (11), we obtain an equation connecting 6 out of 8 
values of / at the vertices of an elementary cube, of the type (xi)' 



f: 



13 



fl2 fl2_ 



/23 fl3_ 



f 



13 



h 



0. 



(14) 



Thus, there naturally appears the decomposition of the cubic lattice into tetrahedra (which 
support the linear problem) and octahedra (which support the nonlinear equations). 

Each equations in (13) can be interpreted as a conservation law and allows for an 
introduction of a potential. For instance, due to the first equation in (13) one can introduce 
a function p on the vertices of the cubic lattice, such that 



92 - Pi, 



P3 ~ Pi, 
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which brings the second equation in (13) into the form 

(P23 - P13)(P3 - Pi) = (P23 - Pl2)(P2 ~ Pi), (15) 

which is equivalent to (x3)- Analogously, the second equation allows us to introduce a 
vertex function q such that 

92 , 93 
a = — , o = — , 
q q 

which brings the first equation to the form 

912-913 /l 1\ 

923 • (16) 



91 v 93 92 

which is nothing but (xs)- Finally, setting 

T\ T\2T T13T 

q = — =^ a = , = , 

T T\T2 T1T3 

we can rewrite equation (16) as 

(7W) (™-?™*) =0 , 

which yields the Hirota equation (xi) 

T2T13 — T3T12 = CT1T23- (17) 

In order to obtain equation (X2)) one has to start with the linear problems 

ip 2 - ip = u(ip! -if)), ip 3 - tp = v(ip! - tp). (18) 

These linear problems are not essentially different from (11): on one hand, they are gauge 
equivalent, and on the other hand, they have already appeared in (12), so that on lattices 
of a sufficiently large dimension both types of linear problems coexist. The condition of 
compatibility of (18) leads to the nonlinear equations 

(u 3 - l)(v - 1) = (v 2 - 1)(«- 1), U 3 Vl = V 2 Ul, (19) 

which are interpreted, as before, as a map (u,ui,v,v{) t-> (^3,^2). Variables tp satisfy 
equation of the type (X2): 

{ipl - 1pl 2 ){lp2 ~ ^23)(^3 ~ ^13) = _ 1 / 2Q ^ 
(tpl2-th)(lp23-i>3)(tpl3-1pl) 

Like in the previous example, both equations (19) can be interpreted as conservation laws 
and allow for an introduction of a potential. However, this time these two equations are 
similar and lead both to the same equations. For instance, one can resolve the second 
equation in (19) by introducing a vertex function / according to 

h h 

u = y, V = T> 
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and then the first equation in (19) leads again to (14). 

One can add one more equation which appears via the substitution 

h = — = , 

<?3 T1T23 

which leads to 

(h 12 - l)(h 3 - 1) = h 2 h 13 (l - hj 1 )^ - fi£). (21) 

This equation will be discussed in Section 7 under the name (Y). One can express the 
variables h directly through tp, by composing all intermediate transformations: 

4>2 - V> _ 1 

~ 93 ~ b ~ h - h ~ 1/v - 1 ~ j>i-j> _ x ~ (^2 - ^3)^1 - ^) ' 

It can be checked directly that this substitution turns equation (21) into a product of 
four copies of equation (20), at the original vertex and at its shifts in three coordinate 
directions. 

Thus, all equations (xi)— (xs) an d (21) are related to one another in a rather simple 
manner: 

(15) = ( X5 ) (21) = (Y) 

p h 
I t 

ip — > u, v <— f — > a,b q r 

(20) = (X2) (14) = (x 4 ) (16) = (x 5 ) (17) = (xi) 



3 Multidimensional consistency 

3.1 Consistency for equations of the cube type 

The idea behind the definition of multidimensional consistency given in Section 1 (Defini- 
tions 2, 3) is rather general and can be implemented for various types of discrete systems. 
However, it refers to the properties of equations on the whole multidimensional lattice, 
which is difficult to verify. It would be preferable to have some sufficient conditions for 
multidimensional consistency which refer to some local (finite) piece of the lattice. 

For 3D equations of the cube type, like dBKP equation (7), such a local sufficient 
condition refers to one 4D cube, as follows (see, e.g., [3]). A solution of a non-degenerate 
equation of the cube type which is generic on each 3D sublattice can be defined by the 
initial data consisting of the values of x on all two-dimensional coordinate planes. Indeed, 
these data are obviously independent, and the inductive application of the equation allows 
one to extend the solution from the coordinate planes to the whole of Ij N (compute 
from the known x, x% and Xij), provided one does not encounter contradictions in this 
inductive process. It is sufficient to verify the lack of contradictions within one 4D cube, 
which is done as follows. Initial data within one 4D cube are: 

x, Xi, x^ (1 < i < j < 4). 
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Application of equations on the four cubic faces adjacent to the vertex x yields the values 
x-ijk (l<i<j<k<4:), and then application of equations on the four cubic faces 
adjacent to the vertex X1234 yields four a priori different values for this last field. The 
4D consistency takes place (and implies multidimensional consistency) if these four values 
identically coincide (thus, one has three conditions in terms of 11 initial data to be fulfilled). 
A more detailed discussion can be found in [3], where it was shown that the dBKP equation 
satisfies this criterium. 

3.2 Lack of consistency for the dKP equation on the face centered cubic 
lattice 

The type of initial value problem discussed for the cube type equations in the previous 
section is not applicable to 3D equations of the octahedron type, because of the absence 
of the fields Xijk from the equations (10). Hence, the very notion of the multidimensional 
consistency in this concrete situation has to be modified. 

Here we analyze a possible definition of 4D consistency for the dKP equation suggested 
by its original form (1). This definition turns out to be unsuccessful, so that this section 
is not necessary for the further reading, however we hope that it will clearly demonstrate 
the non-triviality of the problem of finding the suitable notion. 

Equation (1) decomposes into two independent systems, one for the fields x : Z 3 ven — > C 
on the so called even (or black) sublattice 

Z 3 ven = {m = (mi,m 2 ,m 3 ) G Z 3 : mi + m 2 + m 3 = (mod 2)} , 

and another one for the fields x : Z 3 dd — > C on the odd (or white) sublattice 

^odd = { m — (mi, 7712,7713) G Z 3 : mi + m 2 + 7713 = 1 (mod 2)} ; 

it will be enough to consider the half defined on Z 3 ven (say) . The latter is known as the face 
centered cubic (fee) lattice. Its set of vertices is clearly in a one-to-one correspondence to 
Z 3 , but its (Delaunay) cell structure is quite different. Its 2-cells are equilateral triangles, 
while its 3-cells are octahedra and tetrahedra. Equation (1) relates six fields assigned 
to any elementary octahedron of Z 3 ven , centered at m G Z 3 dd . Two octahedra with a 
nonempty intersection either share an edge (those whose centers are neighbors in Z 3 dd ) 
or either share a vertex (those whose centers are diagonal neighbors in Z 3 dd , i.e., are at 
distance 2 from one another). 

One could attempt to define the multidimensional consistency of octahedron type equa- 
tions by imposing them on all Z 3 vcn sublattices in 

^even = { m = ■ • • > m Jv) G % N '■ 771 1 + . . . + = (mod 2)} . 

However, this idea turns out to be invalid. Indeed, one can show the inconsistency of two 
copies of (1), corresponding to the sublattices (123) and (124): 

x{m + e 3 )x(m — 63) — ax(m + ei)x(m — ei) + bx(m + e2)x(m — e 2 ) = 0, (22) 
x(m + e4)x(772. — 64) — cx{m + e\)x{m — e\) + dx(m + e2)x(m — e 2 ) = 0, (23) 
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where a, b, c, d are arbitrary constants. As initial data for these two equations one can take 
the values of x at four two-dimensional planes parallel to the coordinate plane 12: 

x(mi,m 2 ,0, 0), x(mi,m 2 , -1,0), x(m\, m,2, 0, — 1), i(mi,m 2 ,-l,-l). 

These data are free, in the sense that they are not subject to any equation among (22), 
(23), or remaining two equations corresponding to the sublattices 134 and 234, since any of 
these equations contains at least one pair of points differing by 2e3 or by 2e 4 . In particular, 
consider the following points in these two planes: 

(mi,m 2 ) = (±2,0), (1,±1), (0,0), (-1, ±1), (0,±2) with (m 3 ,m 4 ) = (0,0); 
(mx,m 2 ) = (±1,0), (0,±1) with (m 3 ,m 4 ) = (-1,0), (0,-1); 
(mi,m 2 ) = (0,0) with (m 3 , m 4 ) = (— 1, -1). 

These are round points on Fig. 2. Now equations (22), (23) allow us to compute x at the 
triangular points 

(0,0,1,-1), (0,0,-1,1), (±1,0,1,0), (0,±1,1,0), (±1,0,0,1), (0,±1,0,1), 

and then two different values at the square point (0,0, 1, 1). It can be checked that, for 
any choice of the non-vanishing coefficients a, b, c, d, these two values of x(0, 0, 1, 1) do not 
coincide identically (as functions of the initial data). 



Figure 2. Consistency check for equations (22), (23) leads to a contradic- 
tion at the uppermost vertex. 

Note that two neighboring octahedra in different three-dimensional lattices 1? with 
two common coordinate directions (as those underlying equations (22), (23)) do not share 
triangular faces; rather, they share two pairs of antipodal points, or an 2D equatorial 
square-formed section spanned by these points. 
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3.3 Consistency of the dKP equation on a 4D cube 

Here we prove a local statement about the 4D consistency for the Z 3 version (6) of the 
dKP equation, which deals with the octahedra contained within one elementary 4D cube 
of Z 4 . Each of the eight 3D cubic faces of the 4D cube contains such an octahedron, so 
that we have the eight octahedra 

[j/cTTl] — {xj,X/ s ,X rn ^Xji t ,Xj rn ,Xl tm \ and Ti[jklTl\ — %iki £im, %ijki ^ijrai %ikm\ 

(where {i, j, k, m} = {1, 2, 3, 4}). Neither of the equations on these octahedra involves the 
fields x and £1234, so that only 14 out of 16 vertices are involved. It turns out that one 
can take 9 of them as independent initial data, for instance, 

Xl, x 2 , X 3 , £4, X U , 2 2 4, ^34, £134, 2:234, 

and to find the remaining 5 fields using only three of the equations, namely those for the 
octahedra [124], [134], [234]: 

21224 — 214x2 + 22421 = 0, 

21324 - 21423 + 23421 = 0, (24) 

22324 - 224^3 + £34^2 = 0, 

and for their shifted copies T 3 [124], T 2 [134], Ti[234]: 

£123^34 - £134^23 + £234^13 = 0, 

£123^24 - £124^23 + £234^12 = 0, (25) 
£123^14 - £124^13 + £134^12 = 0. 

We use equations (24) to determine 212, 213, 223. Then, we use the first equation in (25) 
to determine 2123. Finally, we have two alternative answers for 2124 which come from 
the last two equations in (25). We now show that these two values coincide as functions 
of the initial data, so that the 4D consistency takes place. The following computations 
might make an impression of a repeated application of a certain skillful trick, but we will 
show later that this trick is actually a key structural feature common for all consistent 
octahedron type equations, namely the so called tripodal form of these equations. 

First of all, we show that the values of 212, 213, 223 determined from (24), together 
with the initial data x\, 22, 23, automatically satisfy the dKP equation on the octahedron 
[123]. To this end, we rewrite (24) as 



212 


£14 


+ 


£24 


= 


2l2 2 


2124 




2 2 2 4 




£13 


£14 


+ 


£34 


= 


2i2 3 


£l£4 




2 3 2 4 




£23 


£24 


+ 


£34 


= 


£2^3 


£2£4 




2 3 2 4 





(26) 



An obvious linear combination of these equations immediately leads to 

£12 _ £13 £23 _ 

£l£ 2 £l£3 £2£3 
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which is equivalent to 

X12X3 - x 13 x 2 + X23X1 = 0. (27) 

Second, we show that the values of X12, X13, £23 determined from (24), together with 
X14, X2A, X34, automatically satisfy an equation which literally coincides with dKP. For 
this aim, we rewrite equations (24) in another equivalent form: 

X12 X2_ X\ _ Q 

X14X24 X4X24 X4X14 

^-^ + ^ = 0, (28) 
X14X34 X4X34 X4X14 

X23 X 3 X 2 _ Q 



X24X34 X4X34 X4X24 

A suitable linear combination of these equation leads to 

X12 X13 + X 2 3 _ Q 



X14X24 X14X34 X24X34 
which is equivalent to 

X12X34 - X13X24 + X23X14 = 0. (29) 

Finally, we show the 4D consistency. The first equation in (25) used to determine X123 
is equivalent to 

^34 Xi34_ + ^34_ = Q 

X13X23 X13X123 X23X123 

The last two equations in (25) used to determine X124 are equivalent to 

X24 X124 , x 2 34 _ n 



X12X23 X12X123 X12X123 f oi\ 
X14 _ X124 X134 _ Q \ aL ) 



X12X13 X12X123 X13X123 

They give the same value of X124 if and only if 

Xl4 X 2 4 X134 X 234 _ ^ 



X12X13 X12X23 X13X123 X23X123 

Combining the latter equation with (30), we see that the 4D consistency condition is 
reduced to 

Xl4 X 2 4 + X34 _ Q 



X12X13 X12X23 Xi 3 X 2 3 

which is equivalent to the already proven equation (29). 
3.4 Cell structure of the lattice Q(A^) 

Considerations of Section 3.3 give a hint towards a valid formulation of the notion of 
multidimensional consistency for the octahedron type equations. In particular, equation 
(29) for the even shifts calls for an interpretation as an equation on a further octahedron 
present in the 4D lattice. Such interpretation is enabled by an embedding of the lattice 
Q(As) into the root lattice Q(A^), given in (9). We start with a short description of the 
Delaunay cell structure of the four-dimensional lattice Q{A±) [10, 29]. We do not go into 
detail and give an elementary description appropriate for our purposes. For each N there 
are N sorts of iV-cells of Q(A N ) denoted by P(k, N), k = 1, . . . , N. 
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Two sorts of 2- cells: 

P(l,2): black triangles {xi, Xj, x k }, 
P(2,2): white triangles {xij,Xik,Xjk}; 

Three sorts of 3-cells: 

P(l, 3): black tetrahedra {x^, Xj,x k , xg}, with all four facets being black triangles, 

P(2,3): octahedra [ijki] = {xij,Xik,Xi£,Xj kl Xji,X}-e}, the eight triangular facets 
of each octahedron being bi-colored, consult Fig. 3; 

P(3,3): white tetrahedra {x^^, Xy£, x^, Xjke}, with all four facets being white 
triangles; 

Four sorts of 4-cells: 

P(l, 4): black 4-simplices {x^, Xj,Xk, X£, x m }, with all five facets being black tetra- 
hedra, 

P(2,4): 4-ambo-simplices, in the terminology of [10], with the ten vertices Xij, 
i,j £ {0,1,..., 4}; each such polytope has five octahedral facets (i) = 
[jk£m] and five black tetrahedral facets Ti{xj, x k , X£, x m }; 

P(3, 4): 4-ambo-simplices with the ten vertices Xijk, i, j, G {0,1,2, 3, 4}; each of 
these polytopes has five octahedral facets Ti[jk£m] and five white tetra- 
hedral facets; 

P(4, 4): white 4-simplices {x ijke , x ijkm , x ijem , x ikim , x jMm }, with all five facets be- 
ing white tetrahedra. 



jk 




im 



Figure 3. Vertex enumeration of the octahedron (?) = [jk£m], where 
i G {0, 1, 2, 3, 4} and {j, fc, £, m} = {0, 1, 2, 3, 4} \ {i}. Opposite vertices 
carry complementary pairs of indices. Facets are bi-colored; there are 
four white triangles like {jk,j£, k£} missing one of the indices (m in this 
case), and four black triangles like {jk, j£, jm} = Tj{k,£,m} sharing 
one common index (j in this case). 
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The combinatorial description of higher dimensional cells in the root lattices Q(An) 
with bigger N is given analogously. 

The affine Weyl group is generated by permutations of indices and translations. It 
permutes the cells of the same sort. 

We finish this section with a suggestive description of the admittedly somewhat com- 
plicated combinatorics of a 4-ambo-simplex P(2,4) whose five octahedral faces (0), (1), 
(2), (3), (4) carry a quintuple of consistent octahedron type equations. Remarkably, this 
description not only has a combinatorial meaning, but also has a direct relation to the 
multidimensional consistency of equation {x%)-> explained in [7, p. 285]. Consider a map 
x : Q{An) — > MP n satisfying the following condition: the image of any white triangle 
} is a collinear triple of points. Such maps were introduced in the three- 
dimensional situation by Schief [ ] under the name "Laplace-Darboux lattices". He also 
observed the relation of their four- dimensional consistency to the Desargues theorem (pri- 
vate communication; see [2, 6]). These maps (under the name "Desargues maps") are 
studied in detail in the recent work by Doliwa [12, 13]. A connection of the Desargues 
theorem to equation (%4) appeared in [1]. 

It is easy to realize that the image of an octahedron [ijkm] is a complete quadri- 
lateral, as on Fig. 4. It contains four lines which are images of the white triangular 
faces {x{j,Xik,Xjk}, and four triangles which are images of the black triangular faces 
Ti{xj, Xk, x m }. An analytic description of the complete quadrilateral is given, according 
to the classical Menelaus theorem, by equation (x?)- 



01 




Figure 4. Geometric interpretation of equation (x2)'- a Menelaus con- 
figuration, or a complete quadrilateral, is an image of the octahedron 
[0123]; four lines correspond to the white triangles {jk, j£, M}, while 
four triangles correspond to the black triangles {jfc, j£, jm}. 

The image of the 4-ambo-simplex P{2, 4) is then a configuration like the one on Fig. 5. 
One can recognize here five complete quadrilaterals, corresponding to the octahedra (i) = 
[jk£m], as well as the images of the five black tetrahedra Ti{xj,Xk,xi,x m }. The ten lines 
are the images of the ten white triangular faces of the octahedra. The configuration on 
Fig. 5 illustrates one of the most important incidence theorems of the classical projective 
geometry - the Desargues theorem. Its five-fold symmetry is not obvious from the first 
glance (and from the original formulation of the theorem), but it was well known to the 
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classics, see, e.g., [19]. 




Figure 5. The image of a cell P(2, 4) - Desargues configuration IO3. 
Ten lines correspond to the the white triangles {jk, j£, k£}. One can 
clearly recognize the images of the five black tetrahedra {ij, ik, i£, im} = 
Ti{j, k, I, m}, and five complete quadrilaterals, which are images of the 
octahedral faces (i) = [jkim], each one has six vertices missing one of 
the indices (i in this case). 

Remark 4. The image of the 4-ambo-simplex P(3,4) turns out to be less interesting. It is 
a configuration IO254 in the projective plane, i.e., five lines in general position with their 
ten pairwise intersection points, see Fig. 6. This figure does not support any non-trivial 
incidence theorem. 

3.5 Consistent triples of octahedron type equations 

A local definition of the 4D consistency of octahedron type equations should deal with the 
4-cells of Q(A/l) possessing octahedral faces, i.e., with 4-ambo-simplices. We consider in 
more detail an ambo-simplex whose tetrahedral faces are black. Its five octahedral faces 

(0), (1), (2), (3), (4), 

are characterized by the property that any two of them share a triangular face. Two 
4-ambo-simplices are said to be adjacent if they share a common octahedral face. For 
instance, for any fixed i = 0, 1, . . . , 4, the octahedra 

ToTr^O), TxTf^l), T 2 Tr\2), Tgir 1 ^, T^r 1 ^) 

serve as the faces of a 4-ambo-simplex adjacent to the previous one, with the common 
octahedron (i). 

In the rest of this section, we will work with the representation of the lattice Q{A±) 
as Z 4 , which corresponds to the representation of the octahedron type equations in the 
form (8). This representation appears through a projection from Z 5 to Z 4 along the 0-th 
coordinate axis (that is, by forgetting the index no). It is useful to remember that the 
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Figure 6. The image of a cell P(3,4) - a configuration IO254. Five 
lines correspond to the the white tetrahedra {ijk,ij£,ik£, jk£}. One 
can also recognize five complete quadrilaterals, which are images of the 
octahedral faces T{ (i) , each one has six vertices sharing one of the indices 

(i in this case). 



shift operators Tj in the Z 4 representation of Q{A^) stand for the operators T^Tq 1 in the 
standard representation (9) of Q(A±). As already mentioned, we will consider triples of 
octahedron type equations corresponding to three 3D sublattices of Z 4 , say the sublattices 
(124), (134), (234): 

F(xi,X 2 ,X4,Xi2,Xi4,X 2 4) = 0, 

G(xt,X3, X 4 , Xl 3 ,Xu, £34) = 0, (32) 
H(X 2 ,X 3 , X 4 , X 2 3,X 2 i, £34) = 0. 

An elementary combinatorial structure supporting these equations is a triple of octahedra 

(3), (2), (1), (33) 

which are just three out of five octahedral faces of a 4-ambo-simplex (in particular, any 
two octahedra of a triple share a triangular face). We will say that another triple of 
octahedra of the same coordinate directions is adjacent to the triple (33) if they are faces 
of an adjacent 4-ambo-simplex, and the common octahedral face of the 4-ambo-simplices 
belongs to neither of the triples. Thus, for two adjacent triples there is a further (seventh) 
octahedron sharing a face with any octahedron of the both triples. We will say that the 
both triples of octahedra are adjacent along the seventh one. For instance, the triple 

T 3 (3), T 2 <2>, Ti(l), (34) 

is adjacent to triple (33) along the octahedron (0), and the triple 

T 3 T^(3), T 2 T 4 - 1 (2), T^il), (35) 

is adjacent to triple (33) along the octahedron (4). We are now in a position to formulate 
a local definition of 4D consistency of the octahedron type equations. 



17 



Definition 4. A triple of the octahedron type equations (32) is called J^D consistent if they 
can be imposed (admit generic solutions) on any two adjacent triples of octahedra. 

Thus, Definition 4 requires to consider equations (32) on the triple of octahedra (33) 
and on its two adjacent triples, (34) and (35). 

To verify the first requirement of Definition 4, one can consider the same set of 9 initial 
data as used in Section 3.3 for the dKP equation, namely 

Xl, X 2 , 2 3 , 2 4 , 2 M , 2 2 4, £34, 2134, 2234- (36) 

Equations for the three octahedra (33) determine the fields 212, 213, 223: 

£12 = /(21, 2 2 , 24, 214, 224), 

£13 = #(21,23,24,214,234), (37) 

£23 = h(x 2 ,X 3 ,X4,X24,X U ), 

then equation for the octahedron T 3 {3) determines the field 2123: 

£123 = f(xi3, X23, X34, %iu, £234), (38) 

and finally two equations for the octahedra T 2 (2) and Ti(l) deliver two a priori different 
values for 2124: 

£124 = £l(xi2, £23, £24, 2123, £234) 

you J 

= h(x 12 , 213, ^14, ^123, £l34)- 

The first requirement in Definition 4 is that these two values of 2124 identically coincide 
as functions of the 9 initial data (36). 

One can proceed similarly to verify the second requirement of Definition 4: consider 
the set of 9 initial data 

2i, 2 2 , 2 3 , 2 4 , 214, £24, £34, £13,-4, £23,-4- (40) 

One starts with equations (37) for the three octahedra (33), solves equation for the octa- 
hedron T 3 T^ 1 (3) to determine the field 2123,-4: 

£123,-4 = 7(£3, £13, £23, £13,-4, £23,-4), (41) 

and finally has two equations for the octahedra T 2 T^ 1 (2) and TiT 4 ~ 1 (l) which give two a 
priori different values for 212,-4: 

£12,-4 = g{x 2 , 212 , 2 2 3, £23,-4, £123,-4) , . 

u \ (42) 

= h[X\, 212, 213, £13 -4, £123,-4) ■ 

The second requirement in Definition 4 is that the two values of 212,-4 coincide as functions 
of the 9 initial data (40). 

For the formulation of the following statements we use the following convention: the 
derivatives are denoted by the lower indices: 

df dh 

J l — 7^ , • • ■ , «34 



o , ■ ■ • , ' "01 o 
C72l OX34 

(since we will not need higher order derivatives, the multiple indices should cause no mis- 
understanding). 
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Proposition 3. If the triple of the octahedron type equations (37) is J^D consistent, then 
the following holds: 

(fl h /A (h fu fa \ 

rk 5i g 3 gA < 2, rk g 4 g u g u < 2. (43) 
\0 h 2 h 3 h^J \/i 4 /l 2 4 ^34/ 

Proof. The first claim follows from the consistency of equations for the triples (33), (34), 
while the second follows from the consistency of equations for the triples (33), (35). Each 
of the matrices in (43) contains the derivatives with respect to the fields not belonging 
to the octahedron along which the corresponding triple is adjacent. Both statements are 
verified similarly, therefore we only prove the first one. We differentiate condition g = h 
in (39) with respect to the initial data: 

dxi ■ 512/1 + 5123/1351 = ^12/1 + hi 3 gi + ^123/1351 > 

9 X2 '■ 512/2 + 523 ^2 + 5123/23^2 = ^12/2 + ^123/23^-2, 

dxa ■ 523/13 + 5123 (/l3ff3 + /23/13) = /11353 + ^123(/l353 + /23/13), 

d X4 '■ 512/4 + 523^4 + 5l23(/l354 + /23/14) = /112/4 + h 13 g 4 + Zll23(/l354 + f23h 4 )- 

These equations say that the row vector 

(512 — hl2, —hi 3 + (5123 — hi2 3 )fl 3 , 523 + (5123 - ^123)/23) 

belongs to the left kernel of the matrix 

h h / 4 \ 

51 53 54 J , 
h 2 h 3 hi/ 

so that the latter matrix has rank < 2. □ 

Proposition 4. If the triple of the octahedron type equations (37) is consistent, then 
some octahedron type equations are automatically fulfilled on the sublattices (4) = (123) 
and (0) = (1234): 

K(xi,x 2 , X 3 , X12, X13, X23) = (44) 

and 

L(xi2,Xl 3 , X 23 , Z14, £24, ^34) = 0. (45) 

Proof. The fact that 

(fl /2 /A 

rk 5i 53 54 < 2 
\0 h 2 h 3 hj 

can be reformulated as follows: if we consider (37) as a system for the unknowns x\, X2,x 3 ,x 4 
and solve the first two equations of this system for x\, X4, then the substitution of the re- 
sult into the last equation will cancel out X2,x 3 identically. This is precisely statement 
(45). Relation (44) follows similarly. Thus, consistency of a triple of equations on two 
adjacent triples of octahedra yields that the vertices of the connecting octahedron also 
satisfy certain equation. □ 
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Thus, the 4D consistency of a triple of the octahedron type equations yields that 
actually all coordinate directions are on the same footing: each of the five root lattices 
Q(A 3 ) contained in Q(A^) carries its own equation of the octahedron type. 

Remark 5. It is not clear how to extend the initial data (36) or (40) in order to get a set of 
independent initial data for the whole 4D lattice Z 4 . A possible formulation applicable to 
the whole of Z 4 is the following: the initial data consist of the values of x on the coordinate 
2D planes (14) for all i = 1,2, 3 (they intersect along the coordinate axis 4). Within two 
neighboring 4D cubes, one has the following set of 11 independent initial data: 



Xl, X 2 , X 3 , X4, X 14 , X 24 , ^34, X44, X144, X 2 44, X344 



(46) 



(the last four data coming from the second 4D cube). One can determine X12, £13, x 23 
from equations (37), and then X124, X134, X234 from equations (37) shifted in the 4-th 
coordinate direction: X124 = / (xi4, X24, X44, £144, £244) = T±(f), etc. The 4D consistency 
conditions appear from the comparison of the three values of X123 from (38) that must 
coincide identically as functions of the initial data: 

Z123 = f{g,h,x 3 4,T4(g),T4(h)) = g(f, h, x 2 4, T 4 (f), T 4 (/j)) = 5, xu, ?4(/), T 4 (5)). 

(47) 

We will show that this notion of 4D consistency quickly leads to the same necessary 
conditions formulated in Proposition 3. Denote the exterior functions f,g,h in (47), as 
well as their derivatives, by the bar, so that, for instance, 

f = T 3 (f) = f(g,h,x 34 ,T4(g),T 4 (h)), /1 = T 3 (/i). 

The following equations are obtained by differentiating (47): 



dx 3 
dxA 



d x 
d Xl 

d. 
a 



£244 
I 

Z344 



fm =9ih = h 2 fi + h 3 gi, 
M2 = 51/2 + 93^2 = h 2 f 2 , 
fi93 + hh 3 = g 3 h 3 = h 3 g 3 , 
f\9A + hhi = 51/4 + 53^4 = ^2/4 + h 3 g A , 

/l474(54) + /24T 4 (/l 4 ) = ffl4?4(/4) + 534?4(/l4) = W^^) + h^T^g^) 

71474(514) = 5i 4 r 4 (/i4) = h 2 4T4(fu) + huT4(gi4), 

/24?4(/l24) = 5l4?4(/24) + 934^4^24) = ^24?4(/24), 
/u^Cfl^) + hiT4(h 3 4) = 534^4(^34) = 7l34?4(534)- 



(48) 



(49) 



Now equations (48) say that the row vectors 

{-9i, /i, /2-53) and (gi-h 2 , -h 3 , g 3 ) 
belong to the left kernel of the matrix 

7i h / 4 N 
91 53 54 

h 2 h 3 /i 4 , 



20 



Analogously, equations (49) say that the row vectors 

(-914, fu, fu ~ 934:) and (514 - h 2i , -hu, 534) 
belong to the left kernel of the (T4-shifted) matrix 

/4 /l4 /24 
94 914 534 
v /l4 /l24 ^-34, 

Thus, both matrices have rank < 2. 

4 Tripodal forms of the octahedron type equations 

We now derive further analytic consequences of the necessary conditions of 4D consistency 
formulated in Proposition 3, which can be formulated as 

/i53^2 + /25i^3 = 0, /253/14 = f 493^2 + fights, (50) 

and 

/l4534^24 + /24514/134 = 0, /24534^4 = f 4934^24 + f2494h 3 4, (51) 

respectively. 

Proposition 5. If the triple of octahedron type equations (37) is J^D consistent, then it 
can be cast into the form 

a(xi,X4,xu) - b(x 2 , X 4 , X24) = p(xi2, X14, X24), 

c(x 3 , x 4 , x 34 ) - a(x\, x 4 , x M ) = q(x 13 , x i4 , x 34 ), (52) 
b(x 2 , x 4 , X24) - c(x 3 , X 4 , X34) = r(x 2 3, X 2 4, X34), 

and simultaneously into the form 

A(x 1 ,x 4: , x i4 ) - B(x 2 , x 4 , X24) = P(xi, x 2 , X12), 

C(x 3 , x 4 , x 34 ) - A(x±,x 4 , xi 4 ) = Q(x 1 ,x 3 , xi 3 ), (53) 

B(x 2 ,XA, X24) ~ C(x 3 , X 4 , X 34 ) = -R(x 2 , X 3 , X 2 3)- 

Proof. We prove that the general solution of equations (50) is of the form 

/ = 0(a(xi, x 4 , X14) - 6(x 2 , x 4 , X24), X14, X24), 
g = ip(c(x 3 , X4, x 34 ) - a(xi,Xi, x i4 ), xi 4 , x 34 ), 
h = x(K x 2, X4, X24) - c(x 3 , x 4 , x 34 ), x 24 , x 34 ), 

which is obviously equivalent to (52). First equation in (50) implies 

/1 5i/53(xi,x 3 ,x 4 ,xi4,a;34) a(xi, x 4 , x M ) 
— (xi,x 2 ,x 4 ,xi4,X24) = - t 77 , r = r 

72 n<2/Mx2,X3,X4,X24,X34) P(x 2 ,X 4 ,X24) 
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(it is sufficient to choose some fixed values of X3 and X34 in the middle expression). Now, 

ag-3 fihz 

= -j— = -7(^3,^4,^34), 

91 h 2 

where 7 denotes the common value of both ratios. Clearly, a, (3, 7 are defined up to a 
common factor, possibly depending on X4. We arrive at 

Pfi + a/2 = 0, ag 3 + 751 = 0, jh 2 + f3h 3 = 0. 

Choosing some a = a(x\, X4, xu), b = b(x 2 , X4, X24), c = c(x 3 , X4, X34) so that a = 01, 
P = b 2 , 7 = C3, we come to the conclusion that the functions f,g,h can be represented as 

/ = 0(a(xi, x 4 , x M ) - b(x 2 , X4, X24), ^4, x u , X24), 
gr = -0( c (x 3 , x 4 , x 34 ) - a(xi, x 4 , x i4 ), X 4 , X14, X34), 
/i = x(6(x 2 , x 4 , x 24 ) - c(x 3 , X 4 , X34), X 4 , X 2 4, ^34)- 

Substituting this into the second equation in (50), one finds: 

hw + %-°< <54) 

where prime denotes the derivatives with respect to the first arguments of the functions. 
Differentiating this equation with respect to x\ and x 2 yields: 



£(x 4 ), 



where 5 is the common value of all expressions. Now, integration yields: 

-i = S(a - b) + A(x 4 , xu) - fi(x4, X24), = 5(c - a) + 1/(2:4, x 34 ) - A(x 4 , x u ), 

-i = S(b - c) + yu(x 4 , X24) - ^(^4, 2:34) 
X 

(the form of the integration constants follows from (54)). At this point we use the remaining 
freedom in the definition of the functions a,b,c, which can be changed to 

a = k(x/i)a + £(x4, xu), b = k{xi)b + m(x4, X24), c = /c(x4)c + n(x4, X34), 

with arbitrary k,£,m,n. Denoting 

4>(a — b, x 4 , X14, X24) = 4>(a — b, X4, X14, X24), 
ip(c - a, x 4 , X14, x 34 ) = V>(c - a, x 4 , x u , x 34 ), 
X(b - C, X 4 , X 2 4, ^34) = x(b - c, X 4 , X 24 , X34), 

it is easy to see that 

(j) = k<p, (j)^ = 4>4 + (k'(a — b) + £4 — m4)4>', 
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so that 

fa 4>a k' £4 m 4 

and analogously for ip, \- The choice k'/k = 5, £4 = kX, 7714 = kfj,, = kv4 leads to 

4>4 = i>4 = Xa = 0, 

and this gives the desired representation. The solution of the second pair of equations (51) 
is of the same structure and leads to representation (53). □ 
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Figure 7. Tripodal forms (52) of equations (37) on the octahedra 
(3),(2),(1), and the tripodal form (55) of equation (45) on the octa- 
hedron (0) sum up to zero 

We call equations (52) and (53) tripodal forms of equations (37). For instance, the 
three terms in the first equation in (52) correspond to the three (triangular) legs (1,4,14), 
(2,4,14), and (12,14,24) of the tripod with the head (4,14,24), see Fig. 7. Adding the three 
tripodal forms (52) of equations on the octahedra (3), (2), (1) leads to the equation 

P(xi2,xi4, X24) + q{xi3,xu, x 34 ) + r(x 2 3, X 2 4, X34) = 0, (55) 

which is nothing but the tripodal form of equation (45) on the octahedron (0) , with the 
head (14,24,34). 

Similarly, adding the three tripodal forms (53) of equations on the octahedra (3), (2), (1) 
leads to the equation 

P(x 1 ,x 2 , x 12 ) + Q(x 1 ,x 3 , xiz) + R(x 2 , x 3 , x 23 ) = 0, (56) 
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Figure 8. Tripodal forms (53) of equations (37) on the octahedra 
(3),(2),(1), and the tripodal form (56) of equation (44) on the octa- 
hedron (4) sum up to zero 



which is the tripodal form of equation (44) on the octahedron (4), with the head (1,2,3). 
Thus, Proposition 5 subsumes (and is actually much stronger than) Proposition 4. 

Moreover, there exist further tripodal representations of equations (37) and (44), (45), 
due to the symmetry of all coordinates. Each equation under consideration admits eight 
tripodal representations, since each face of the octahedron can be chosen as the head of the 
tripod. In total we have 40 such representations. In order to put them in a unified form, 
we will use the realization of the Q(A±) lattice as a hyperplane (9) in Z 5 , and will denote 
the leg functions in the tripodal form by the letter a with three superscripts which are the 
indices of the three arguments, omitting the arguments themselves. Our convention will 
be that this notation is symmetric with respect to the first and the third argument (the 
base of the leg) , while the second argument (the spike of the leg) enters the equation only 
once. In this notation, formulas of Proposition 5 can be written so: for the quadruple of 
octahedra (i), (j), (k), (m), the tripodal forms 



(i>: 


jn,jm,mn 


^kn,km,mn 


„ jn,jk,kn 

a J J , 




<J> = 


^kn,km,mn 


^in,im,mn _ 


kn,ik,in 
a , 


(57) 


(k): 


^in,im,mn 


„ jn,jm,mn 


. a in,ij,jn 





sum up to the tripodal form 

(m) : a jn ' jkM + a kn ' lk ' in + a in ^ jn = 0. 
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Figure 9. Left: enumeration of the vertices of an octahedron: (1,2,3) is a 
face, and opposite vertices carry complementary indices which sum up 
to 7. Right: the tripodal form with the head (i,j,k). 



The following theorem summarizes the results of the present section. 

Theorem 6. A quintuple of the octahedron type equations is J^D consistent if and only 
if, in any quadruple of octahedra such that any two of them share a triangular face, the 
tripodal forms of any three equations sum up to the tripodal form of the fourth one. 

Proof. The necessity follows from Proposition 5, the sufficiency is proved by literally the 
same argument as used in Section 3.3 for the proof of the 4D consistency of the dKP 
equation. □ 



5 Tripodal equations 
5.1 Definitions and notation 

We have shown that each equation of a 4D consistent system can be written in 8 ways in a 
tripodal form. In the present section, we temporarily forget about consistency and analyze 
just this property of a single octahedron type equation. The enumeration of the vertices 
of an octahedron and of the corresponding variables will be different from the previously 
used: it will be convenient to enumerate them just from 1 to 6 as on Fig. 9. 
Our task here will be to describe all equations 

$(x 1 ,X2,X 3 ,X 4: ,X 5 ,X6) = 0, (58) 

which are locally equivalent to equations of the form 

a(xi,x 7 _ k ,Xj) + b(xj,x 7 -i,x k ) + c(x k ,x 7 -j,Xi) = 

for each of the eight triples (i,j,k), corresponding to the faces of an octahedron. Such 
equations will be called tripodal. 
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As usual, it is supposed that equation (58) is irreducible, in particular, neither of 
the partial derivatives <I>j vanishes identically. The answer is determined up to point 
transformations 

Xi = Xi(xi), X'^O. (59) 

Functions a, b, c may be different for different faces. It will be convenient to denote them 
by superscripts pointing to the arguments, for instance, a 123 = a 123 (xi, X2, x 3 ). Here the 
second index is distinguished, and it is supposed that the notation is symmetric with 
respect to the first and the third index: a 123 = a 321 . This means that o 123 (xi, X2, £3) = 
a 321 (x3, X2, x\) (but of course it is not supposed that a (xi, X2, £3) = a (x\, X2, x 3 )). 
Opposite vertices of an octahedron will be denoted by the same letter in the lower and the 
upper cases, e.g., I = 7 — i. Thus, the tripodal form of equation with the head (ijk) takes 
the form 

a iKj + a jlk + a kJi = 

see Fig. 9. The set of all eight tripodal forms looks like this: 



head: 


tripodal form: 




(1,2,3) 


a 142 + a 263 


+ a 351 


= 0, 


(1,2,4) 


a 132 + a 264 


+ a 451 


= 0, 


(1,3,5) 


a 123 + a 365 


+ a 541 


= 0, 


(1,4,5) 


a 124 + a 465 


+ a 531 


= 0, 


(2,3,6) 


a 213 + a 356 


+ a 642 


= 0, 


(2,4,6) 


a 214 + a 456 


+ a 632 


= 0, 


(3,5,6) 


fl 315 + a 546 


+ a 623 


= 0, 


(4,5,6) 


a 415 + a 536 


+ a 624 


= 0. 



Clearly, this structure is rather restrictive. It turns out to be possible to find a complete 
list of tripodal equations. It is given and discussed in the following subsection, while the 
rest of the section is devoted to the proof that the list is exhaustive, indeed. 

5.2 Classification of tripodal equations 

Theorem 7. All tripodal equations (58), up to the point transformations (59) and trans- 



positions i -<->• I, (i,j) (J, I), are given by the following list: 

x\x G + X2X5 + £3^4 = 0, (Ti) 

(x 1 - X4){x 2 - Xq)(x 3 - X 5 ) + (X4 - X 2 )(x 6 - X 3 )(x 5 - X\) = 0, (T 2 ) 

(x\ - x 2 )x A + (x 2 - x 3 )x 6 + (x 3 - xi)x 5 = 0, (T 3 ) 

xix 6 = (x 2 + x 3 )" y (x 4: + x 5 ), (T 4 ) 

XlX 6 = x 2 + x 3 + x 4 + x 5 , (T 5 ) 

XlX 2 X 3 X 4 = X 5 + XQ, (T 6 ) 

XI + X2 + X 3 + X4 + X5 + XQ = 0. (T 7 ) 



It is easy to see that the first three equations are equivalent to (xi), (X3)> (X2), upon 
re-naming (X12, £13, £23) — > (^4, X5, x%). The proof of the theorem will be given in the 
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Eq. 


legs 


case 




xyz 


T 
1 


(Ta) 


\y + zJ 


III 


(T 3 ) 


(x + z)y, log I ), log(x + y) 


III 


(T 4 ) 


xy, y(x + z) 1 ' wfx + zW 7 , logfx + y), w 


I (7 / 1), III (7 = 1) 


(T 5 ) 


y, (x + y)z 


II 


(T 6 ) 


a;y, y, y + log(x + z), log(x + y) 


I 


(T 7 ) 


y 


III 



Figure 10. Legs store 



following subsections. Here we will just check that all equations of the list are tripodal, 
indeed, that is, admit all eight representations (60). 

The results of this check are summarized in Table 10, which lists the functions a(x, y, z) 
acting as the legs a Jjfc = a(xi, Xj,Xk) for equations of the list. Recall that the variables x, z 
are on the equal footing, while y plays a distinguished role. In particular, the functions a 
always depend on y but is sometimes independent on x or z. Functions a are defined up to 
point transformations (59). Besides, they can be multiplied by an arbitrary constant, one 
can flip x, z and add arbitrary combinations of the type n(x) + v{z) (as long as one such 
function is considered, and not all three legs of the tripodal form). At this stage, we are 
not concerned with adjusting all these transformations for consistent equations, therefore 
the table contains arbitrary (possibly simple) representatives for the leg functions. This 
arbitrariness notwithstanding, we will see in Section 6 that this table is very useful for 
putting separate tripodal equations into consistent quintuples. The third column of the 
table contains the class according to the subdivision at the end of Section 5.3. 

Equation (Ti). The tripodal form with the head (1, 2, 3) is 



:i,2,3) : 



x,i 



+ 



X6 



+ 



a; 5 



0. 



X\X2 X2X3 X3X1 



other ones are obtained by reflections Xi «-)■ xj and (xi,Xj) f> (xj, xi), which generate the 
symmetry group of an octahedron and leave the equation invariant. 

Equation (T2). The multiplicative tripodal form (1,2,3) can be seen directly from 
the equation. All other ones follow again by reflections Xi 4-¥ xj and (xi,Xj) (xj,xi), 
although the invariance of the equation under these reflections is less obvious than in the 
case of equation Ti. 

Equation (T3) is already written in the tripodal form (123). Arranging the terms in a 
different order leads to the tripodal form for the opposite face: 

(4, 5, 6) : xi(x4 - x 5 ) + x 3 (x 5 - x 6 ) + x 2 (x 6 - x 4 ) = 0. 
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Thus, the equation is invariant under the involution 

P : (xi,x 2 ,x 3 ) o (x 6 ,x 5 ,x 4 ). 
Further, one finds the multiplicative tripodal form 

(1,2,4): (g^) (*4-* 5 ) =L 

\X 3 - X 2 / (X6 - X 4 ) 

All other tripodal forms are obtained with the help of the involution P and the cyclic 
permutation 

Z : X\ — >■ X2 — >■ X 3 — > Xl, Xq — > Xq — > X 4 — )• X6- 

Equation (T 4 ) is invariant under the transformations 

xi^x 6 ; (x 2 ,x 5 ) o (x 3 ,x 4 ); (xi, x 2 , x 3 , 7) o (x^ 1/7 , x 5 , x 4 , -I/7), 
which allows one to obtain all the tripodal forms from the following two: 

(1.2.3) : ^- Xe ( X2 + X3 y + ^ = 0i 

Xl X\ 

(1.2.4) : (logxi +7log(x 2 + x 3 )) +logx 6 - log(x 4 + x 5 ) = 0. 
Equation (T5). In this case the symmetry group is generated by the transformations 

Xi <H> xr, (^2,^5)^(^3,^4), 
which allows one to obtain all the tripodal forms from one of them, for instance, from 

X 4 + Xo X 3 + X5 

1,2,3 : -x e + 3 5 = 0. 

Xl Xl 

It is worth mentioning that, although here the first and the third terms are similar to the 
legs encountered for equations (T3) and (T 4 ) with 7 = 1, this is a different type of legs, 
since the role of the distinguished variable y is different. 
Equation (Tq). Here there are transformations 

(xi,x 6 ) O (x 2 , x 5 ); x 3 O x 4 , 

which make it sufficient to give the following three tripodal forms: 

(1, 2, 3) : xix 4 X2 = 0, 

(1, 3, 5) : log(xix 2 x 3 ) - log(x 5 + x 6 ) + log x 4 = 0, 

(3,5,6) : (logx 3 + logxi) + (logx 4 - log(x 5 + x 6 )) + logx 2 = 0. 

Equation (T7). This case is trivial. 



28 



5.3 Reduction to functions of two variables 

Let us prove some useful relations which follow from the (local) equivalence of the two 
tripodal forms 

(*, j, K) : a ikj + a jIK + a KJi = 0, (i, j, k) : a iKj + a jIk + a kJi = 0, (61) 

which correspond to two faces sharing a common edge ij. We will use the subscripts i to 
denote the partial derivatives d Xi - 

Proposition 8. The following identities hold true: 




Since each one contains five variables only, they are satisfied identically and not by virtue 
of equations (61). 

Proof. The numerators do not vanish identically because of the local solvability of equa- 
tions with respect to each variable. We shall prove only the first identities in each pair, 
since the second ones follow by the permutation (i, I) (j, J). Note also that all identities 
are invariant under the flip k <-> K. Solving equations (61) for xj, we obtain an equation 
of the following form: 

f( Xj ,x K , + a KJi ) = g(x 3 ,x k , + a kJi ), 

and by differentiation there follow equations of the form 

fK + fa$ Ji = gfa$i, f(af j + 4 Ji ) = g'(af j + af i ), f >a« Ji = g 'a k /\ 

where the prime denotes the derivatives of the functions /, g with respect to their last 
arguments. Now (62) follows immediately from the second and the third equations, while 
from the first and the third we derive 

f KJi 

IK dj iKj KJi l i ikj i KJi 

-JT = ~JUi a K ~ a K =<t> => <Pja k = 4> k aj , 
J aj 

which yields 

( ( a J Jt \ JKj _ KJi) ikj _ , KJi\2 JKjf 1 \ 

\\WJj k J k ~ WJ k 

Dividing by a 1 ^ ^ a %kj '(aj''') 2 and differentiating with respect to Xj, we arrive at (63). □ 

Equations (62) will play the main role in the subsequent analysis. Equations (63) are 
needed only for the proof of the following property, which, in its turn, will only be used in 
the proof of Proposition 18. Nevertheless, it is convenient to do the job right now, since 
all necessary formulas are at hand. 
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Proposition 9. If at least one of the functions a//^ , a^j 1 , a^ k k , a k j l does not vanish 
identically, then 



ikj ikj iKj iKj 

a ik a kj _ a iK a Kj 

ft n n A rv>K 



(64) 



Proof. Differentiating (62) with respect to x and x , we obtain equations 



n KJi - u iK \ n kJi ' u kj I jlK J _ "Kj I j/fc 

Along with (63) they build a pair of linear homogeneous systems, one for (1 / 'aj Jl ) k , 
{l/a k j J% )k and another for (1/ 'a?/ )k, {^/^/ k )k- Determinants of these systems coincide. 
The statement of the proposition says that if at least one of the systems has a non-trivial 
solution, then their common determinant vanishes. □ 

The following representation is a direct consequence of identities (62). 

Proposition 10. The functions a* fcj and a lK ^ are of the form 

a ikj = a ij b k + p ik + p kj^ a iKj = a ij b K + p iK + p Kj^ b k b K ^ Q ^ 

Proof. Differentiating equations (62), we find: 

n ik 3 n kJi _ JKj KM ikj j Ik _ iKj UK 

a ij a J — a ij a J ) a ij a I — a ij a I ■ l DD J 

Setting the variables xi,xj,xk to arbitrary constants, we find: 

af/ = a tj f3 ik = a ij j jk af/ = a ij b k , 

and an integration leads to formula (65). Here the function a 13 is defined up to addition 
of the terms /i 1 + v 3 and multiplication by a constant, which can be taken into account 
by a re-definition of b and p. It is easy to see from (66) that this function can be chosen 
the same in the both formulas, and that the factors b can be taken non-vanishing without 
restriction of generality. □ 

Remark 6. Taking into account all the tripodal forms (60), we come to the conclusion that 
each function a lmn has a representation of the type (65), so that the form of the equation 
is already found up to functions of two variables. However it should be understood that 
the notation in formulas (65) refer to some fixed pair of the tripodal forms (61). If we 
would like to apply them to all possible index sets simultaneously then we would be forced 
to use more complicated enumeration for the functions b and p. However, there will be no 
need to do this, since we will only perform a pairwise comparison of the tripodal forms. 

Proposition 11. If a?- ^ then either the both functions b k and b K are different from 
constants or the both are constants, the tripodal forms (61) being either 

a ij b k + ( — + ( — = & a l3 b K + ^- + — = Q (67) 
b K b K b k b k 
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or 

Jk | kj\ i (41 , ifjN | / Ji , iK 



^ (a*' + p** + p^') + (e^ 7 + p fcj ) + (c Ji + p lk ) = 0, 

respective/?/. 

Proof. We still did not exhaust the content of identities (66). If ^ then they are 
reduced to 

b k a kJi _ b K a KJi = Q) h k a 3lk _ h K a jIK = g 5 (6g) 

which give, upon integration, 

akJi = C -j^+ dkl > a3lk =^ + dJk > ( 7 °) 
a KM = ~ + d™, = ^ + d? K . (71) 

Thus, the tripodal forms (61) can be put as 

( a vft* + p * +p^ ) + + + + = 0, 

(atfft* + p « + p^) + + ^ fc ) + + ^) = 0. 
It is easy to see that they are equivalent if and only if 

b K (p ik + p kj + d jK + d Ki ) = b k (p iK + p Kj + <i J ' fe + d fci ), 

from which it follows b% (pf +p k k j ) = b k k (pf+p^ j ). Now if b\ = 0, then pf + p k k j ^ 0, 
since otherwise the first tripodal form would not contain x k , but then also b^ = 0, and, 
setting b k = b K = 1, we find formulas (68). 
If, on the contrary, b k 7^ 0, then 

p ik + p k i = & fc (p* + p>) +q i + gj, p iK + p^ = 6^(p* + pJ) + ? + 

It is not difficult to see that these terms can be assumed to vanish (this can be achieved 
by a re-definition of other terms). But then 

b cr 

and these terms can be again absorbed by the functions c? 1 , c 7 *. As a result, the tripodal 
forms (65) can be put as (67). □ 

The further analysis does not require for complicated computations, but rather for 
a detailed bookkeeping of different cases. Taking into account the possibilities listed in 
Proposition 11, one can take care of this by separating the following classes of equations: 

I. a^A 7^ for at least one triple i, k,j; 

II. a k J- = for all i, k,j, but cif? 7^ for at least one triple; 

III. d k ? = for all i, k,j. 

Clearly, this exhausts all logical possibilities. A complete description of these cases will be 
sufficient for a proof of Theorem 7. 
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5.4 Class I 

Setting, without restriction of generality, b k = Xk, b K = xk, we put the tripodal forms 
(67) as 

a ij x k + — + — = 0, a ij x K + — + — = 0. (72) 

XK XK Xk X k 

Now we have to compare them with the other tripodal forms. Since functions a^ Ik , a kJi 
have to be of the form (65), as well, we immediately find that 

cP 1 = (9b 1 + d 1 , c Ji = cV + d J . 

Assume first that Cj ^ or cf l ^ 0. For definiteness, let cPj ^ 0. Compare (72) with 
the tripodal form k): a 3 + a? %k + a kJI = 0. Apply formulas (70), from which there 
follows: 

j 

a jik = ^_ h i +g ji +q jk^ 

Xk 

while the function a iK i should admit the both representations 

a iK i = a ij x K = + & a lK3 = 

b l b l W 

Since, by assumption, / 0, there follows b\ ^ 0, but then there holds also b\ / 0, and 
the equation is reduced to the form 

r Kj r kj 

-W + aH ' + jr = °- 

Subtracting this from (72), we come to 



c jI c Ji iklI c kJ 
Xk x k b l 



which yields that the equation is reduced by the point transformations (59) to the form 

(Ti). ' . 

Now let cP = cf l = 0, then equation (72) takes the form 

log a lj + log x K - log(c/ + cj) + log x k = 0. 

This case is covered (after an obvious change of enumeration) by the following proposition 
which refers to somewhat more general equations. This is done in order to avoid unneces- 
sary repetitions for the cases II and III, which also lead to equations of this type. In this 
proposition it is not assumed that the equation falls into one of the three classes separated 
above. Actually, the analysis of results performed in Section 5.2 shows that equation (Tg) 
belongs to the class I, equation (T7) belongs to the class III, and equation (T4) belongs to 
the class I if A 7^ 1 and to the class III if A = 1. 

Proposition 12. Let one of the tripodal forms of equation (58) be as follows: 

p IK + p Ki + pl k + q j + g J = Q (?3) 

Then it is reduced to one of the equations (T4), (Tq), or (Tf). 
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Proof. Setting, without restriction of generality, g J = xj, q J = xj, we put the equation as 

p iK + p Ki + pl k + x . +XJ = _ ( 7 4) 

Let us show that if / 0, then = and = 0. 

In order to prove the first of these equations, compare two tripodal forms (61), 

(t, j, K) : a ikj + a jIK + a KJl = 0, (i, j, it) : a^' + a j/fc + a fc Ji = 0, 

assuming that (74) is the first one of them, so that a ifcj = p lk + Xj, a? 1K = p IK and 
qKJi _ pKi _|_ From the first identity (62) we find: 

JKj , kJi iKj 
ik , Ki _ a i + a i Ki _ a iK f7 c\ 

a J a J 

Suppose that p*g\ 7^ 0- Then, differentiating the latter equation with respect to xj, 
we find a k j l = a kJ j, but then the first equations implies also a^j\ % = 0. Therefore, a kJl = 
[ix j + r kl , so the tripodal form k) becomes 

a iKj + a jlk + ^ xj + r ki = Q ^ ^ Q 

Eliminating xj with the help of (74), we obtain the identity 

a iKj + a jlk _ ^ p IK + p Ki + p ik + +r ki = Q ^ p IK = Q _ 

The second equation is shown analogously, by comparing (74) with the tripodal form 
(I,j,K) (it is enough to use the symmetry i K, k -H- I). 

Thus, we have shown that (74) takes one of the following two forms: 

p 1 + p Ki + p k + xj + x j = or p IK + p ik + Xj + xj = 0. 

But it is clear that the first case is included into the second one, upon the re-labeling 
i—^k—^I—^K—^i. Therefore we can assume that in equation (74) there holds p K% = 0. 
Note that this makes the equation invariant under the re-labeling i -f* k, I •<-)• K, which 
does not affect the tripodal form (i,j,k). Now, returning to identity (75), we find: 

iKj , kJi 

ik _ a i + a i JKj _ iKj _ n 



Pi = JUi a iK = a 



and, using the symmetry just pointed out, 

Jlk , kJi 

ik _ a k + a k ilk _ jlk _ n 

Pk ~ hji a Ik ~ a jk - U - 



As a result, we put the tripodal forms (i, j, K) and k) as 

p iK + p ik + Xj . + XJ = 0i + a ji + a fcJ< 
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and, solving for xj, we arrive at the identity 



IK Kj u 



p IK + Xj = ^{x k ,x h a^ + a^) - P_]^ = a _i^ = * 



Pi 



l K J 1 a 1 



where the last equality is obtained by setting X j = const. Upon suitable point transfor- 
mations of xj and Xk, we can put p IK as p IK = P(xj + xk), while = a?/ = rK There 
follows easily 

Kj , jl 
1 = a j +a J 

P' r i ' 

and, differentiating twice with respect to xj or to xk, we find (1/P')" = 0, or 
P(y) = alog(y + /3) + 7 or P(y) = ay + /3. 

Of course, for p %k analogous formulas are found (by comparison with the tripodal form 
(I,j,K) and upon suitable point transformations of Xi,Xk), so p %k = Q(x{ + Xk) with 

Q(y) = Alog(y + /x) + v or Q(y) = Xy + fi. 

Different combinations lead, after further obvious point transformations, to the cases (T4), 
(T 6 ), (T 7 ). □ 

5.5 Class II 

Now let a^jy = for all i,k,j, but 7^ for at least one triple. For this triple we find 
two equivalent tripodal forms (68) 

(a lj + p lk + p kj ) + (c 3 ' 7 + p K i) + (c Ji +p lK ) = 
O (a ij + p iK + p Kj ) + (c 3 ' 7 + p kj ) + (c Ji + p ifc ) = 

which have to be compared with other ones. Here it is convenient to distinguish subcases 
depending on which mixed derivatives of the pairs of functions p lk ,p k i and p lK ,p •> vanish. 

First assume that for at least one pair the mixed derivatives of the both functions do 
not vanish. For the sake of definiteness, let p] k / 0, pA 7^ 0. Note that formula (69) with 
constant b k ,b K yields that 

Applying this implication to the second of the tripodal forms, we find: 

Jib _ n jlk _ kj , n kJi _ iKj _ n Ji _ %K _ n . 

a jlk — U ' a jk — Pkj T u => a Ji — a iK — U =^ c Ji — PiK — u > 

n kJi _ n kJi _ ik / n , jlk _ iKj _ n JI _ Kj _ n 

a kJi — u ; a ki — Pik 7^ u =^ a j/ — a Kj — U ^ S'/ ~~ ~~ ' 

and as a consequence the equation takes the form 

a ij + + + ? / + ^ + ^ = 0; 
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and after the change q n — > x n we come to equation 

p ij + pjk + p ki + XJ + x j + XR = q ( 76 ) 

If the previous subcase does not take place, then we assume first that one of the pairs 
pik^pkj Qr ptK^pKj con t a ins exactly one function with vanishing mixed derivatives. For 
definiteness, let p %k k ^ and p fe j = 0. Then, as before, = p 1 ^- = 0, and we come to an 
equation of the form 

a ij + p ik + q 1 + c Ji + p iK = 0, 
or, upon a point change of variables, to equation 

pij + pik +p iJ +p iK +Xi = Q ^ 77 j 

Finally, if all mixed derivatives Pi^PkpPiK^PK'j vanish, then our equation belongs to 
the special type already dealt with in Proposition 12: 

a ij + q k + d 1 + q K + c Ji = 0; 

we have seen that such equations cannot belong to the case II. 

As a result, case II is reduced to equations of two special types. They are dealt with 
in the following two statements which result in equation (T5) as the only possible one. 

Proposition 13. There exist no equations of the class II with one of the tripodal forms 
as in (76). 

Proof. We compare (76) with the tripodal form 

(/, J, K) : a IkJ + a JiK + a KjI = 0. 
Solving for X}~, we obtain an identity 

f(xi, xj, a J%K + a KjI ) = g{xi,Xj,xi + xj + x K ), 

whence 

log af tK - log af jI = log gi - log 5^- = h(x i: Xj,xi + xj + x K ). 

Differentiating with respect to xj, xj, and xk, and using the symmetry of the equation 
under permutations of i,j, k, we find: 



a KjI 
a jl 


JiK 
a Ji 


JiK 
a iK 


a KjI 

a Kj 


Kjl 

V 


n JiK 
a i 


n JiK 
u i 


a KjI 


n IkJ 
a kJ 


a KjI 

Kj 


Kjl 


n IkJ 
a Ik 


n Ik.J 
a k 




Kjl 
a 3 


n IkJ 
a k 


n JiK 
a iK 


41' 


IkJ 
a kJ 


n JiK 

a Ji _ 


„JiK 
u i 


„IkJ 
a k 


n IkJ 
a k 


n JiK 
u i 



x 1 , 
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where A , fJ. J , and v K denote the common values of the corresponding expressions. Since, 
by assumption, the equation belongs to the class II, we have a J j^_ = 0, and, differentiating 
the last equality in the first equation with respect to xj, we find: 

„ JiK „ JiK 
J K = _ a iK a Ji = q 

and analogously v K \ = X 1 fi J = 0; since there holds additionally A + fi J + v = 0, all 
three functions actually vanish. This means that the tripodal form (UK) is reduced, upon 
a point transformation of Xi,Xj,Xk, to the form analogous to (76): 

a IJ + a JK + a KI + Xi + xj + x k = 0. 

Resolving once again the both tripodal forms for x^, we find an identity 

a IJ + a JK + a KI + Xi + Xj = g(x i} x^xi + xj + xk), 

whence 

a \J + af 1 = a 1 / + a J / = a J K K + af 7 . 
These equations are easy to solve, and we arrive at the tripodal form 

\{xi + xj + xr-) 2 + /i(:rr + xj + x^) + z/ + Xi + Xj + = 0. 

The case A = leads to equation (T7) which does not belong to the class II, while for 
A ^ 0, as easily shown, this tripodal form cannot be equivalent to equation (76). □ 

Proposition 14. Equations of the class II with one of the tripodal forms as in (77) can 
be reduced to (T5). 

Proof. Comparing with the tripodal form 

(/, J, K) : a IU + a JiK + a KjI = 0, 

we find: 



iJ „JiK 1 „IkJ JlK „JiK j_ „KjI 

(78) 



Pj _ a J + a J Pk _ a K J a K 

rrAk „IkJ ' ij h 

Pk a k PS a j 



There follows aj 1 ^ = 0. Further, assume that ctji 7^ 0. Then, differentiating the first 
equation with respect to xj and Xj, we find a^ J = ajj J = 0, so that the tripodal form is 

a kJ + a JiK + a K ^ = 0. 

Solving for x/, we find: 

p V + p ik + p U + p iK = f( XKjXjja kJ + a MK^ (79) 

and applying djdk leads to = /j, where the prime stands for the derivative of / with 
respect to the third argument. But then p|j = 0. 
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Exchanging the roles of j and k, we prove that the implications hold: 

JiK / n . IkJ _ AkJ _ ij _ n . JiK i n . Kjl _ Kjl _ ik _ n 
a Ji T 1 U => a Ik — a IJ — Pij — U > a iK T U => a jl — a KI — Pik — U - 

If ajj K / and afj^ / simultaneously, then pV. = = 0, and our equation reduces to 

P> + p fe + J + p iK + xi = 0, 

but this is, up to notation, a particular case of the previous special equation (76). Let, for 
example, aj] K ^ and aj^ = 0. Then (79) is fulfilled, and, applying djdx, we come to 
= f K (recall that a j 1 ^ = 0), but then also p]^ = 0, and we come again to an equation 
of the type (76): 

p* + p ik + p iJ + p K + x/ = 0. 

There remains the only possibility cijl K = af^ = 0. In this case there follows from 
(78) that 



pT\ = (pJ_ 

ik j I ij 

,PkJi \Pj , 



0. 



and, using the symmetry of our equation under k O K (that is, comparing (77) with the 
tripodal form {IJk) instead of (UK)), we find also: 



%K I I ij 

J>K Ji \Pj , 



0. 



All these equations yield: 



pij + p ik _|_ p iJ + p iK _ +q k + q J + q Kyi + g i ^ 

moreover rj ^ const, since otherwise we would come to equation (T7) which does not 
belong to the class II. Up to changes of variables, equation (77) is reduced to 

Xj + Xk + Xj + X K = XlXi + s\ 

The tripodal form (I, J, K) in this case is a IkJ + a 1 + a K ^ = 0, and, solving for Xj, we 
find: 

IkJ 1 „i\ . 1 i _ a i 



xiXi + s l - Xk - xj - x K = f(xK,xi,a +a l ) =>■ xj + s 



n ikj- 

a k 



Differentiating the last equation with respect to Xi and xi, we obtain: a % u = 0, but then 
s-j = 0, and an additional linear change of variables leads to equation (T5). □ 
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5.6 Class III 

In this case it turns out to be convenient to change notation and to re-write the tripodal 
forms (60) as 

faces: equations: 

(1,2,3), (4,5,6) p 14 + p A2 + p 26 + p Q3 + p 35 + p 51 = 0, 

(1, 2, 4), (3, 5, 6) q 13 + q 32 + q 2G + g 64 + q A5 + q 51 = 0, (80) 

(1, 3, 5), (2, 4, 6) r 63 + r 32 + r 21 + r 14 + r 45 + r 56 = 0, 

(1,4, 5), (2, 3, 6) s 64 + s 42 + s 21 + s 13 + s 35 + s 56 = 

(through a re-ordering of terms each of these equations gives two tripodal forms, corre- 
sponding to opposite faces). 

Proposition 15. Each term in equations (80), as a function f(x,y) of its arguments, is 
of the following type: 

f = a{x)b(y) + c(x) + d(y) or f = plog(a(x) + b(y)) + c(x) + d(y). (81) 

Proof. First we show that 

f*= r{ ?l\(\ +'vW> f y= ( v^~( s + u ^- ^ 

(3{x)+5(y) n(y) + x(x) 

Due to the symmetry of formulas (80) it is enough to do this for / = p 14 , x = x%, y = x±. 
We use the first identity (62), with (i,j,k) = (1,2,3) and (4,5,6): 

aj 32 + a 451 a} i2 + af l a 465 + a 12i a 415 + af 4 



This yields: 



,451 „351 ' „124 ^624 



q\ 3 + qf p\*+pf sf + sf pf+pf 

qf + qf pf+pf sf + 4 2 Pf+Pf 



(83) 



and, setting all variables except for x\,x^ to constants, we come to (82). 
Further, computing the mixed derivatives, we find: 

a(x)S'(y) X(yW(x) 
Ixy (f3(x) + 5(y)) 2 (My) + x(x))2- 

If at least one of the functions a, 5', A, x' vanishes, then / = c(x) + d(y) . Otherwise, taking 
logarithms and applying d x d y , we come to 

P>(x)6>(y) tUy)*Hx) 



(I3(x) + 5(y)) 2 (p(y) + x(x)) 2 - 

If /3' = fjf = 0, then equations (82) are easily integrated and lead to / = a(x)b(y) + c(x) + 
d(y). Finally, if f}' and fj do not vanish, then the two latter equations yield a//3' = X/ fx' = 
const, and an integration finishes the proof. □ 
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The proposition just proven is not easy to use directly, since it is not yet known how 
are the representations (81) corresponding to different terms of the tripodal forms (80) 
related to one another. We start with filtering away the case when these forms contain 
too many terms of the kind a(x) + b{y). 

Proposition 16. If each of the forms (80) contains terms with vanishing mixed deriva- 
tives, then one of the forms contains at least four such terms, and it is a particular case 
of the equation (73). 

Proof. By differentiation we derive from (83): 

«li = O p\i = ^ qH = 0, 

that is, the presence of one term with vanishing mixed derivatives in one of the forms 
yields the presence of such terms in further two forms. This chain can be continued to a 
closed cycle: 

p\\ = => q\t = =► sit = => Ptl = => qE = => s 2 2 \ = => p\\ = 0, 

not containing r lJ . Of course, all forms are on the same footing, so that starting, e.g., from 
r 14 , we would get the cycle 

r\\ = => pf 2 = =► q 2 2 ! = 0=>rg = 0=>^ = 0=> gf} = =► r\\ = 0, 

not containing s 13 . Therefore, if the terms with vanishing mixed derivatives are present 
in each of the forms, then one of the forms accumulates at least four such terms. For 
instance, if both cycles above take place then p\\ = pg| = pjj = P35 = 0, and the first 
equation in (80) becomes p 4 + p 26 + p 3 + p 51 =0. □ 

To compare the general terms (81) we will use the following lemma. Note that the 
function F in this lemma defines a general solution of the Liouville equation (logF) xy = 
-2F. 

Lemma 17. All solutions of the functional equation 

g x {x,z)h y (y,z) 



(g(x,z) - h(y,z)) 2 

are given by the formulas 



F{x,y)^0 



5{Z) + e( Z ), h= +«(*), F — 



y j(z)-a(x) w ' J(z)-f3(y) w ' (a(x) - /3(y)) 2 ' 

Proof. Integration with respect to x leads to 



h y (y,z) 
g(x,z) - h{y,z) 



Mv, z ) ~ J F(x,y)dx. 
Setting y = yo = const, we come to 



- j(z) - a(x), 



g(x,z) - e{z) 
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and since h y 7^ 0, we can choose yo so that the numerator and the denominator do not 
vanish identically. This yields the expression for g. Further, we can use the invariance of 
the equation under non-degenerate linear-fractional transformations 

p{z)g + q{z) p(z)h + q(z) 

9^~r~\ — , — t~\i h ^ ( \u 1 — rr> f^f 

r{z)g + s{z) r{z)h + s(z) 

to bring g to the form g = a(x), a' 7^ 0. Then 

a' (x)h v (y, z) „, . h vz 2h z h z a' 

\ ) v\vi ) - F(x,y) => -j^ + t Vt =0 => 71. — , v\9 ~ °» 



(ct(x) — h(y, z)) 2 h y h — a(x) (h — a(x)) 5 

that is, h = /3(y), and the inverse linear- fractional transformation finishes the proof. □ 

Now we are in a position to analyze the case where at least one of the tripodal forms 
(80) contains no terms with the vanishing mixed derivatives (representable as a sum of 
functions of a single variable). For definiteness, let it be the first form in (80). 

Proposition 18. Ifp+j f- 1 for all i,j, the equation (58) is equivalent to (T%) or to (T3). 

Proof. We apply Proposition 9. We have: a lk i = p %k + p k ^ , a^ IK = pi 1 +p IK , a KJl = 
pKj -\- p Jl ^ a nd, since by assumption a?^ = p 1 ^ 7^ 0, the following equation is fulfilled: 

ik kj iKj iKj 

PikPkj _ a iK a Kj _ A ij 



Then, according to Lemma 17, we have: 



6 k 5 k 



n ik , kj 

P k ^Pk 7 k_ a i 1 k_pj' 

so that the representations (81) for p lk and for p k i are tied in the following sense: 

if p lk = a l b k + c l + d k , then p k i = a?b k + e 5 — d k ; 

if p tk = p\og{a l + b k ) + c 1 + d k , then p k i = —p\og{a^ + b k ) + cP — d k . 
Applying this to the pair pP 1 , and further cyclically, we obtain the representations 

ptk _ a t^k _|_ c i _|_ ^ pik _ pl g( a « _|_ Jjk^ _|_ c i _|_ ^ 

p k i = a j b k + c j - d k , p kj = -plog(a j + b k ) + c j - d k , 

pi 1 = a^b 1 - c j + d 1 , pi 1 = p \og{a j + b 1 ) - c j + d 1 , 

p IK = a K b I + C K _ d I^ <>' - p IK = _ p l Qg ( a K + b I^ + c K _ d I^ 

pKJ = a K b J _ C K + d J ) p KJ = p l Qg ( a K + b J} _ C K + d J^ 

pji _ d i b J _|_ gi _ d J ^ pji _ _pi g(g« _|_ 5^) + c* — d J . 

A comparison of the first and the last terms shows that 

a l = £a l , & = -c 1 + 5, a 1 = ea 1 + 5, c l = -c % + a. 
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Summing up all the forms, we obtain, after point transformations and re-enumeration, the 
following equations: 

xix 4 + x 4 x 2 + x 2 x 6 + x 6 x 3 + x 3 x 5 + ex 5 xi = 5, 
(x 1 + x 4 )(x 2 + x 6 )(x 3 + x 5 ) = a(x 4 + x 2 )(x 6 + x 3 )(x 5 + exi + 5). 

One can be more precise about the constants here. For this aim, we solve these equations 
for one of the variables, say for xq = f(xi, x 2 , x 3 , x 4 , x$), and compare with the second 
tripodal form (80), which yields another expression for 

x 6 = tp(x 2 ,x 4 , q 13 + q 32 + q i5 + q 51 ). 

It follows that (log(/3//5))24 = 0, which can be checked in a straightforward way for a 
given function /. This check shows that the constants are uniquely determined, and we 
come to the equations 

X\X 4 + X4X2 + X 2 Xq + X G X 3 + X 3 X 5 - X5X1 = 0, 
(xi + x 4 )(x 2 + x 6 )(x 3 + x 5 ) = (x 4 + x 2 )(x 6 + x 3 )(x 5 + Xl). 

The first of them is reduced to (T3) after an additional change x 2 — > —x 2 , xq — > —xq, 
while the second one coincides with (T 2 ) after changing the signs of X4,x^,xq. □ 

This finishes the proof of Theorem 7. 

6 Classification of compatible quintuples 

6.1 Separating away non-compatible equations 

Now we have to combine tripodal equations into compatible quintuples, taking into ac- 
count that the legs on the faces shared by different octahedra must coincide, according to 
Proposition 5. Table 10 allows us to exclude some a priori non-compatible combinations. 
For instance, equation of the type (Ti) can only match either an equation of the same 
type or an equation of the type (Tq). Other combinations are impossible, as shown in the 
following proposition. 

Proposition 19. The functions 

xyz, xy, y, (x + z) k y (k ^ 0), (x + y)z, 
y + log(x + z), log(x + y) -log(y + z), log(x + y), 

are pairwise non- equivalent modulo transformations 

a(x,y,z) = ja(f{x),g(y),h(z)) + /j(x) + u(z) 

with non-constant f,g, h, along with the flip x <-> z. Functions (x + z) k y with different k 
are also non-equivalent. 

Proof. All proofs are similar, therefore we consider several pairs of functions as examples. 



41 



• The equality 

yg(y) + 7^og(f(x) + h(z)) + n(x) + v(z) = log — ■ — , 

y + z 

is impossible, as follows immediately by differentiation with respect to x and z. 

• Suppose that 

lf{x)g{y)h{z) + fi(x) + v{z) = (x + y)z, 
then the differentiation gives jf'(x)g'(y)h'(z) = 0. 

• Suppose that 

7 (/(x) + h{z)) k g(y) + + v(z) = (x + z) m y. 

Then g(y) = ay + /3, and further aj(f(x) + h(z)) k = (x + z) m , so /(x) + h(z) = 
const(x + z) m l k . Differentiation with respect to x, z yields m = k. 

The reader is invited to work out the other pairs of functions. □ 

A more careful analysis shows that equations (Ti) and (Te) cannot be matched, al- 
though they have a common leg. More precisely, the following statement holds true. 

Proposition 20. There do not exist compatible triples of equations where one of the 
equations is of the type (T4) with 7 7^ 1, or (T5), or (Tq). 

Proof The proof is based on the fact that each of the equations listed in the proposition 
possesses a unique leg, which does not appear by equations of other types (y(x + z) 7 , 
(x + y)z, and y + log(x + z), respectively). By virtue of Proposition 5, there follows that 
at least one further equation of the compatible triple has to be of the same type. Since 
all coordinate directions Z 4 are on the same footing, we can assume that the compatible 
triple is given by equations (52), the unique leg being the function a(xi, x&, xu) in the first 
and the second equations. Moreover, since this leg is unique, the equations themselves 
are recovered up to transformations (59), in other words, the functions b(x2,X4,X2i) and 
c(x3, X4, X34) can be determined. It turns out that b — c does not depend on £4, so that 
the third equation in (52) is reducible, in contradiction with our standing assumption. In 
the following detailed analysis X{ stands for an arbitrary non-constant function of x%. 
Equation (T4) with 7 7^ 1. Consider equation 

X1X24 = {X2 + + X14) 7 . 

Its tripodal form with the head (4, 14, 24) is 

X11X4 + X14) 7 — X2/X24 = X12/X24, 

and, comparing with the first equation in (52), we find: 

a(xi,X4,xu) = X\(X4 + X14) 7 + \(xu) + n{x4), 

b(x 2 ,X4,X24:) = X2/X24 + (J.(X4) + K^)- 
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This a corresponds to a unique leg y(x + z) 7 in Table 10, therefore the second equation in 
(52) 

-X3/X34 — Xi(Xi + X14) 7 = — X13/X34. 

Therefore c(x%, 24, X34) = X3/X34 + ^(24) + ^(234), but then the terms ^(24) in the third 
equation in (52) cancel. 

Equation (T5). Consider equation 

X2XH = X\ + X4 + X24 + X12 => ^ — = ^r; — . 

A 14 A 14 

We have 

A^i + X4 

o(iCi, 24, X14) = — h A(x M ) + fjtfa), b(x 2 , 24, 224) = X 2 + /u(x 4 ) + ^(224)- 

A 14 

Due to the uniqueness of the leg a the second equation in (52) has to be of the type (T5), 
as well, and to have the tripodal form 

_ Xj + X4 _ X34 + X23 

X14 Xu 

whence c(x%, 24, 234) = X%+ ^(24) + ^(234) , and the third equation in (52) does not contain 
X4. 

Equation (Tg). Consider equation 

X 1 X 2 X l2 X 2 4 = X A + X U => (log(X 4 + Xi4)-logXi)-logX 2 = log(Xi2X24). 
We have 

a(xt, 24,214) = log(X 4 + X14) - log Xi + X(x u ) + fifa), 
6(22,24,224) = logX 2 + ^(24) + ^(x 24 ). 

Also this a is unique, therefore the second equation in (52) has to be of the type (Tq) and 
to have the tripodal form 

logX 3 - (log(A~ 4 + X 14 ) - log Xi) = - log(Xi 3 X 34 ), 

whence 0(23, 24, 234) = logX3 + ^(24) + ^(234), and, as in the previous case, the third 
equation in (52) does not contain 24. □ 



6.2 Completing the classification 

Note that after removing equations (T5) and {Tq) from Table 10 the leg xy becomes unique, 
however the argumentation like in the proof of Proposition 20 is not possible. To handle 
with the remaining cases, it is not enough to use formulas (52) alone, and one has to refer 
to all tripodal forms, (57) and their shifted versions. 

Before we formulate the final result, let us describe more precisely what transformations 
are allowed to bring the equations to the canonical form. First of all, these are autonomous 
point transformations x —> X(x) (the same at all lattice points). However, if only these are 
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allowed, then the answer will contain many arbitrary constant parameters. It turns out that 
all these parameters are inessential and can be killed by non-autonomous transformations, 
which depend on the lattice point. Generally speaking, such transformations result in non- 
autonomous equations, therefore not all them should be allowed but only those special ones 
which render the transformed equation still autonomous. Their existence is related to a 
certain symmetry of the equations. For instance, if the equation is invariant under the one- 
parameter group of translations x — > x+a, then it admits non- autonomous transformations 
x(i,j, k) — > x(i,j, k) + ai + j3j + 7&, with the transformed equation being dependent on the 
parameters a, f3, 7, which can be used to simplify the result. Similar transformations exist 
in all cases under consideration, and we use them to eliminate all the constant parameters 
(the situation is similar for continuous 3D integrable systems). 

Theorem 21. Any compatible quintuple of irreducible shift-invariant octahedron type equa- 
tions on Q{A^) is equivalent, modulo non- autonomous point transformations, to one of the 
following systems ( different indices stand for the shifts in different coordinate directions; 
recall that in the Z 4 realization of the lattice Q{A^) the shift Tq in the coordinate direction 
can be simply omitted). 

Five equations of type (Ti): 



X%jX km Xi k Xj m Xi m Xj k — 0, 

Five equations of type (T2).' 

{xij X; lk )[x k j X km )(Xj m 

(Xik X k j^)(x km •Ejmjfaim Xij) 



0<i<j</c<m<4; 



(xf) 



-1, i,j,k,me {0,1,2,3,4}; ( X 5 2 ) 



Two different quintuples consisting of four equations of type (T3) and one equation of type 
(T 2 ): ' 

(x ik - Xij)x i0 + (xij - x jk )x j0 + (x jk - x ik )x k0 = 0, i,j, k G {1,2,3,4}, 

Ol2 - X 13 )(x 2 3 - X U )(x 2 4 ~ X U ) (X3X2) 



and 



(Xi 3 - X23)(X34 ~ ^24) (X13 - X X2 ) 



( x lk Xij + Xij x jk + x jk x lk = ^ ijj)A . G { 1>2 ,3,4}, 



X i0 Xjo X kQ 

(x 12 - Xi 3 )(x 2 3 - Z34)(£24 - X U ) 



(X4X2) 



1; 



. {x 13 - X 2 3)(x 3 4 - X 2 4){X13 - X 12 ) 

Three equations of type (T4) and two equations of type (T3) 



X 



x 40 \x j0 x m J 

X\3 ~ %12 , Xu - X 23 X 2 3 - X13 



1 



1 



X12 



+ 

^10 %20 X 30 

X\i - X 2 4 X 24: - X34 : X34 - XU 



X23 



+ 



ije {1,2,3}, 
0. 

0. 



(X5X4) 



X13 
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Proof. General scheme. We start with one tripodal equation, replacing variables x^ by 
yet unknown non-constant functions X^ = Xij(xij). Comparing terms in the tripodal 
forms (57), we are able to completely determine the consistent quintuple, up to ten arbi- 
trary functions X12, ■ ■ ■ , X45. To specify the latter, we use the shifted tripodal forms. One 
can show then that all functions Xij are related to one another via some linear-fractional 
transformations. Point transformations allow us to assume that the functions X^ are 
linear-fractional, with coefficients of different functions being connected by certain rela- 
tions. Resolving these relations, we come to a consistent system containing several free 
parameters. Finally, we get rid of them using non-autonomous transformations. 

Equation (Ti) can only be compatible with equations of the same type, as follows from 
Table 10. Equations (57) with (i,j,k,m,n) = (1,2,3,4,0) read: 

(1): X23 = Q, 

-^03-^04 -Xo2^04 X02X03 

< 2 > : !Tx--Y^ + Y^r = ^ ^ 

^03^04 ^01 A 04 ^01^03 
(3) ^24 X\4. ^ X\2 q 

-X"o2-Xo4 XqiXq^ Xq\Xq2 

Indeed, in the first equation the functions Xij can be freely chosen; this defines the first 
term in the second equation. But then in the second term in (2) the dependence on X40 
is known, while the variables X14 and xoi (which were absent from (1)) enter via arbitrary 
new functions. Continuing in this fashion, we express all terms through ten arbitrary 
functions Xij. As a corollary we get 

(4) : X01X23 — A~02^13 + -^03X12 = 0> (0) : ^12-^34 — X13X24 + -^14-^23 = 0. 

Equations (57) are herewith exhausted. To determine the functions X%j we have to consider 
the shifted tripodal forms. We write them as 

r l( i ): Tl f_^__i^ + _i^ = o 

T 2 (2) : 1 



Xq2 


X03 


+ 


Xqa 


X23X2A 


X23X34 


X24X34 


Xq\ 


-^03 


+ 


X04 


-^13^14 


-^13X34 


X14X34 


Xq\ 


Xq2 


+ 


X04 


X12X14 


X 12X24 


X14X24 



T 3 (3): T 3 — ^ + — ^ = 
\A12A14 A12A24 A14A24 

and the comparison of the legs on the left-hand sides leads to relations like 

aX 30 (x ) X 10 (x ) 



^23(^2)^34(^4) X 12 (x2)Xu(x4) 



It is easy to realize that the terms [i + v can be neglected here, and that the functions 
X12, X13, X23 coincide up to constant factors, and the same is true for the functions 
^01,^02,^03- Because of the symmetry of all the indices, we come to the conclusion 
that all ten functions are proportional. Without loss of generality we can 
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and a direct check shows that the equations are consistent for arbitrary values of a^j. All 
these parameters can be set equal to 1 upon use of the non-autonomous transformation 

x(no,ni,n 2 ,n3,n 4 ) = Q™- nj x(n , n 1 , n 2 , n 3 , n 4 ), 

and the answer takes the form of the quintuple (xf). 

Equation (T4) with 7 = 1. We consider this equation in the form 

-^04(^24 ~~ -^34) = -^23(^02 — Xq 3 ). 
Identifying it with the first equation in (57) with (i,j, k, m, n) = (1, 2, 3, 4, 0), we have: 

„02, 24,04 Y Y 03,34,04 Y Y „02,23,03 V (V V \ 

a — A24A04, a — A34Ao 4 , a — A 23 (A 03 — A02J, 

and the last two equations in (57) read: 

(2) : X 34 Xo 4 - a 01 ' 14 ' 04 = a 03 ' 13 ' 01 , (3) : a 01 ' 14 ' 04 - X 24 X 4 = a 01 ' 12 ' 02 . 



Since the leg xy is unique, these equations are of the type (T4), as well. This determines 
them in the following form: 



X34X04 — A*i4 X04 — Xx3 (A"oi — 



03 J 



01 J 



where K, similarly to X, stand for arbitrary non-constant functions of the corresponding 
variables: Y%j = Yij(xij). Summing up, we find the fourth equation 



(4) : X 12 (Y 02 - Y01) + X 13 (X i -Y 03 ) + X 23 (X 03 - X, 



02) 



0, 



which can be of the type (T3) only. But then the corresponding functions X and Y are 
linearly related, and it is easy to see that we can assume that they plainly coincide (re- 
defining X12, Xis, and X15, if necessary). As a result, we come to the following quintuple 
of equations: 



(1) 
(2) 
(3) 
(4) 

(0) 



(A~24 — AT34)A~o 4 = A~23(A~03 — Xq 2 ), 
(X 3 4 — Xi4)Xo 4 = Xi 3 (Xqi — Xqs), 
(X14 — AT24)Ao 4 = A'i2(A'o 2 — A~oi), 
-^12(^02 — Xqi) + Xi3(Aoi — Xq 3 ) + X 23 {Xq 3 
X\A — X 2 4 X24 — X 3 4 X34 — A14 



02 J 



0, 



(85) 



X12 



X23 



X 13 



0. 



Equations (57) are now exhausted, and the remaining functional freedom has to be reduced 
by using their shifted versions: 



Ti(l) : 


Ti 


X23X03 


— ^23X02 


— Xq4 (X24 


~ X34) = 


r 2 <2) : 


T 2 


^13X01 


— Xi 3 Xq 3 


— Xo4(X 3 4 


- X u ) = 


73<3> : 


T 3 


X12X02 


— X\%Xq\ 


— A~04 (X14 


— X24) = 



0). 
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The comparison of legs on the left-hand sides leads to relations like 



aX 23 (x 2 )X 03 (x ) = X 12 (x 2 )X i(x ) + n(x 2 ) + v(x^), 

where the terms n + v account for the freedom in the definition of the legs. There follows 
that the functions Xi 2 ,Xi 3 ,X 23 are the same up to constant factors, while the relations 
between the functions X15, X 2 § and X35 are affine. Further, equation (4) takes the form 

aX m {xQi){X 2A (x 12 ) - X 34 (xi 3 )) + (3X Q/ t (x 02 ) (X 3i (x 23 ) - X u (x 12 )) 
+X 04l (xo3)(Xu(xi3) - X 24 (x 23 )) = 0, 

and the comparison with the previous form leads to the conclusion that X14, X 2 4, and X34 
are affinely related to X12, and Xq^ is affinely related to X15. Finally, for (0) we obtain 
the representation 

X 02 (x 12 ) - X 03 (xi 3 ) | ^03(^23) ~ X i(x 12 ) | ^ ^oi(xi 3 ) - ^02(3:23) _ Q 
X 4(xu) X a(x 2 a) X 4(x 34 ) 

(it is enough to observe that equations (85) are invariant under the flip 4 -H- and the 
inversion of Xq^, X\ 2 , Xi 3 , ^23)1 arid a comparison with the previous representation shows 
that I/X04 is amnely related to X14. As a result, all functions are expressed through one 
of them (say, X\ 2 ), and without loss of generality we can put them as functions of x as 

X 12 = x, X 13 = px, X 23 = qx, X 04 = l/x, 
X i4 = aiX + bi, X 0i = a/x + di, i = 1,2,3. 

Substitution into the shifted tripodal forms allows us to show that b% = di = and leads 
to equations 

J -^- = Ax ij ( J -\, 1,3 = 1,2,3. 

X 4 \x 0j XQiJ 

They are consistent for all parameter values. The non-autonomous change x(i,j, k, m, n) = 
a\a2a 3 A mn x(i, j, k, m, n) allows us to set all parameters to 1, and we come to the consistent 
triple (x|x!)- 

Equation (T3). Now we do not consider equations of the type (T4) anymore, and the 
leg y(x + z) becomes unique for equations of the type (T3). It is not difficult to see that in 
a consistent quintuple containing such an equation at least four equations are of the same 
type. Indeed, either the tripodal form does not contain legs y{x + z), or all three legs are 
of this sort. Therefore, if one starts with such a leg for one equation of the triple (57), 
then all the legs for the three equations and for their sum (which is the tripodal form of 
the fourth equation) will be of this sort. We easily come to the following system: 



(1) 
(2) 
(3) 
(4) 



-^24(^02 — -X04) + ^34(^04 — Xq 3 ) + X 23 (X { 
-^34(^03 — X04) + -^14(^04 ~~ ^ 01 ) ^ 1 

X\a{Xqi — X04) + ^24(^04 — Xq 2 ) + Xi 2 (X { 

-^12(^02 — -Xoi) + Xi 3 (Xq\ — Xq 3 ) + X 23 (X { 



03 


— -^02) 


= 0, 


01 


— Xq 3 ) 


= 0, 


02 


-*0l) 


= 0, 


03 


— Xq 2 ) 


= 0, 



(86) 



^ , {X14 — Xi 2 )(X 2 4 — X 23 )(X 3 4 — Xj 3 ) _ ^ 

(X12 — X 2 i) (X 23 — X 3 i) (X13 — X14) 
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Here (4) follows from the first three equations in an obvious way, while (0) is derived using 
their multiplicative representation 



Xi4 — Xij _ Xq4 — Xoj 
Xij — Xj4 Xoi — Xq4 

In order to determine the functions X, we compare legs in the shifted tripodal forms. The 
first three equations are re-written as 



Ti(l) 

r 2 (2) 



T\ [Xq2{X2A — X23) + Xq 3 (X2 3 — X34) + Xo4(^34 — X24) = 0) , 
Ti{Xq\{X\z — X14) + Xo3(X34 — X13) + Xo4(Xi4 — X34) = 0) , 
23(^01(^14 — X12) + Xq2(X\2 — X24) + Xq4(X24: — X\ij = 0) , 



and the legs comparison leads to equations like 

aXo3(zo) (^23(^2) ~ ^34(^4)) = Xq\{xq) (Xi2(x 2 ) - Xi 4 (x 4 )) + /u(x 2 ) + v(x±). 

Because of the symmetry of the indices 1, 2, 3, 4, we find that the functions X\2, X\ 3 , X23, 
X14, -^24 j -^34 are aflinely related, and the same is true for the functions Xq\, Xo 2 , X$ 3 , 
Xq±. To establish a relation between the two sets of functions, it is enough to compare the 
legs Ti(a 03 ' 04 ' 02 ) of equation (1) and T (a 13 ' 14 ' 12 ) of equation (5): 

, X m {x 2 ) - X m (x4) . , X 12 {x 2 ) - X u {xi) 

log v ( \ y ( \ = A log y < \ Y ( \ +K x 2) + v{x 3 ). 87) 

Xo 4 (x 4 ) - X 03 {x 3 ) X M (x A ) - X 13 (x 3 ) 

Differentiating with respect to X2, we find: 

X^) = XX' 12 (x 2 ) + , 
^25(^2) - ^45(^4) ^12(^2) - Xia(xa) 

whence Xq4 and X14 are related by a linear-fractional transformation. Without loss of 
generality, we can restrict ourselves to two cases: 

1) Xij(x) = aijx + bij, i, j = 0, 1,2,3,4; 

2) Xij(x) = a^x + bij, i,j = 1,2,3,4, X 0i (x) = Ci/x + di, i = 1,2,3,4. 

A rather tiresome analysis of the shifted tripodal forms shows that in the case 1) there 
are two possibilities: etjj = c^ay, bij = b\2 or a^- = 1, bij = f3{ + and in the case 2) 
o-ij = C{Cj, bij = bi2, di = d\. Non-autonomous dilations and translations allow us to 
reduce these cases to the quintuples (X3X2) and (xixz)- 

Equation (T2). Now we can exclude from consideration equations of the type (T3). 
Then equation of the type (T2) can be a member of a consistent quintuple with equations 
of the same type only. As above, the tripodal forms (57) allow us to show, starting with 
one equation, that a consistent quintuple must be of the form 

(X^ — Xik)(Xkj — Xj tm )(Xj m — Xim) _ j 

\X-ih Xfcj^Xkm Aj m )(Aj m Xij) 
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Comparison of the shifted tripodal forms leads to relations like (87) with all possible per- 
mutations of indices, whence all X^j are related by linear-fractional transformations. When 
we encountered such a situation above, we used the method of undetermined coefficients, 
but, fortunately, now one can dispense with it. Indeed, the same arguments as in Lemma 
17 allow us to deduce from (87) that the linear-fractional relation between X02 and X12 is 
the same as between A04 and X14 or between X03 and X13 (with A = 1 and fi + v = 0). 
This is true for other index sets, as well, so that 

X ik = MijX jk , i ^ j ^ k / i, 

where My are some linear- fractional function. There follows MijMjk = Mik, and, setting 
Li = Moi, Lq = id, we obtain My = LiLj 1 . But then Lj 1 Xn t = Lj 1 Xjk, in particular, 
L^ 1 Xoi = Lj 1 Xoj, and on the other hand, Xy = L^Xqj = LjXoi, whence L^Lj = LjLi. 
Thus, Xij(x) = LiLjX(x), where all Lj pairwise commute, and X(x) is an arbitrary 
function. Setting without loss of generality X(x) = x, and performing a non-autonomous 
transformation x(no, n\, ri2, rt3, 774) = L™ 1 L^ 2 L^ 3 L™* {x(uq, ni, 712, n^, Ui)), the situation is 
reduced to the case when all Li = id. 

Equation (T7). A simple analysis of the tripodal forms (57) and their shifted versions 
shows that consistent systems of this type are reduced to linear ones with constant coef- 
ficients, which is of course not very interesting. Note however that in this case it makes 
sense to relax the condition of the shift invariance, which might lead to auxiliary linear 
problems for hypothetical 4D equations. This interesting possibility remains outside of the 
scope of the present paper. □ 



7 Remarks on the Y-system 

We are aware of one integrable equation on Z 3 which formally belongs to the octahedron 
type equations and is not in our list. The so called Y-system (of type A jy) 

NI , N (1 — h(m + e2))(l — him — e 2 )) ,^ rS 

him + e^him-ex) = ^ (Y) 

(1 — h i (m + e3))(l - h i (m-e 3 J) 

plays a prominent role in the theory of integrable systems [40, 5, 11, 15, 17, 21, 26, 37, 38]. 
The Y-system is closely related to equations (xi) _ (xs)- However, its properties are rather 
different from those of (xi) _ (xs)- It h as been already shown in Section 2 that on Z 3 this 
equation can be derived, through suitable changes of variables, from both (xi) an d (X2)- 
We discuss here in more detail its relation to equation {x2j- 

Start with a 3D Desargues (or Laplace-Darboux) map as defined in Section 3.4, i.e., 
with a map x : Q(A^) — > IP" such that the image of any white triangle is a collinear 
triple of points. In the representation of Q(As) as Z 3 ven , this can be seen as a 2D Q-net 
(indexed ny mi, 1712) and its iterated Laplace transforms (indexed by 777,3), see [33, 12]; 
one can find definitions of a Q-net and of its Laplace transforms in [7]. Thus, the vertices 
of any white tetrahedron are collinear, and therefore one can introduce their (real) cross- 
ratio, called h (or 1 — h, or else 1 — h , depending on the ordering of the arguments of 
the cross-ratio). Any black tetrahedron has six white neighbors, adjacent to it along its 
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six edges. A straightforward computation shows that the variables h for these six white 
tetrahedra are related by equation (Y) which basically says that the product of six cross- 
ratios is equal to 1. To obtain it, one has to consider four octahedra adjacent to the 
black tetrahedron along its black faces. They carry four multi-ratio equations (X2) coming 
from the corresponding four complete quadrilaterals. These equations, being multiplied, 
give the Y-system (the product of four multi-ratios can be re-arranged as a product of six 
cross-ratios). To sum up: variables h of equation (Y) live on white tetrahedra, while the 
equation itself is assigned to black tetrahedra of Q(A^). 

One can try to make this picture multidimensional, by starting with a multidimensional 
Desargues map x : Q{A^) — > M.P n . Then any 3D sublattice x : Q{A^) — > M.P n will support 
its own copy of equation (Y), with the variables h which show up in this copy and in no 
other: they are defined on 3-cells, and each 3-cell belongs to only one 3D sublattice. The 
only relation between these co-existing copies of equation (Y) will be the condition that 
the product of h's for all tetrahedral faces of a white simplex is equal to 1. 

Thus, although equation (Y) is closely related to multidimensionally consistent octa- 
hedron type equations, there is no reason for itself to be multidimensionally consistent. 
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